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I hereby certi^' that this correspondence is being facsimile transmitted to the U.S. Patent and 
Trademark Office, PCT Legal Office, on the date shown below. 
Attention: Leonard Smith, PCT Legal Examiner 

Anthony Smith, Attorney-Advisor Facsimile: 703-308-6459 

Pages - 5 Pages 

c'%Jc^ /^J^jLm.^^^ Date: August 26, 2002 

Diane Schwaiger u 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 
In re application of 

Ali Jalali, Pierre Leray, and Dominique Art Unit: Not Assigned 
Lacroix 

Application No.: 09/581,272 Examiner: Not Assigned 

Date Received by PCT/PTO: June 7, 2000 Docket No.: 136.147 

Customer No. 023907 

For: Method of Calculating the Fast Fourier Transform and the Inverse Fast 
Fourier Transform 

RENEWED PETITION UNDER 37 CFR 1 .47(a^ 

Assistant Commissioner for Patents 
Box PCT 

Washington, D.C. 20231 

Attention: Mr. Leonard Smith, PCT Legal Examiner 
Mr. Anthony Smith, Attorney-Advisor 

Dear Messrs. Smith and Smith; 

The Decision on Petition Under 37 CFR 1.47(a) mailed June 25, 2002 states that item 
(2), factual proof that the missing joint inventor refuses to execute the application or cannot be 
reached after diligent effort, has not been satisfied. Otherwise, all other items have been 
satisfied. 

In the Decision) the examiners state that "a copy of the certified mail return receipt 
postcard has not been included as an exhibit to the present petition." Applicants' representative 
submits that the original certified mail return receipt postcard showing that the Febraary 1, 2001 
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February 1, 2001 letter to Mr. Jalali was received was attached to ±e Declaration and 
Statement of Facts by Christian Hamon in Re 37 CFR 1.47(a) and Section 409.03(d) MPEP 
that was filed with the Renewed Petition Under 3 CFR 1. 47(a) on April 11 , 2002. A 
photocopy of the certified mail return receipt postcard is enclosed herewith. Also enclosed is a 
copy of our postcard showing that the original certified return receipt postcard was received by 
the PCT/PTO on April 23, 2002. 

As Mr. Hamon stated in the Declaration and Incitement of Facts by Christian Hamon in 
Re 37 CFR 1.47(a) and Section 409.03(d) MPEP, filed on April 11, 2002. he could not 
identify the signature of the person who signed the return receipt postcard. 

Forthcoming Evidence 

As documentary evidence that every effort has been made lo obtan Mr. Jalali' s 
signature on the application papers an additional Declaration and Statement of Facts with 
copies of the following documents will be filed as soon as the executed Declaration is received 
from Mr. Hamon: 

1) French-language letter addressed to Mr. Jalali dated September 28, 2000. an 
English translation of the letter, a Declaration and an Assignmem that was enclosed with the 
letter, and the address label; 

2) French-language letter addressed to Mr. Jalali dated February 1, 2001, an 
English translation of the letter, a Declaration and an Assignment that was enclosed with the 
letter, the address label, and the return receipt postcard signed by a person whose signature is 
not readable; and 

3) French-language letter addressed to Mr. Jalali dated June 26, 2001, an English 
translation of the letter, a Declaration, an Assignment, and a complete copy of the application 
Serial No, 09/581,272 and Preliminary Amendment that were enclosed with the letter. 

The above letters and enclosures are additional documentary evidence that every effort 
has been made to obtain Mr. Jalali's signature on the application papers. None of the letters 
were returned as undeliverable and Mr. Hamon knows of no other means with which to contact 
Mr. Jalali. 
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Re<]|uest for One-month Eartension of Time 
A one-month extension of time from August 25 , 2002 to September 25 , 2002 is 
requested in which to file the additional Renewed Petition Under 37 CFR 1.47(a) and 
Declaration and Statement of Facts by Christian Hamon. Mr, Hamon is on holiday until 
August 26. 2002 and cannot be contacted before that date, I was advised by my French 
associate on August 21, 2002 that Mr. Hamon is expected to sign and return the additional 
Declaration during the week of August 26, 2002, and that it will be sent to me by courier 
shortly thereafter. Mr. Hamon's Declaration will then be filed in the U.S. Patent Office. 

The Commissioner is hereby authorized to charge payment of any fees required to 
Deposit Account No. 14-1080. 

Respectfully submitted. 



Date: August 26. 2002 

Nilles & NiUes, S.C. 
U.S. Bank Center, Suite 2000 
777 East Wisconsin Avenue 
MUwaukee. WI 53202 
Telephone: 414-276-0977 
Facsimile: 414-276-0982 

(b G:\Data\CLIENT\136\147\FTO-RMponse-8.doc 




James £. Nilles 
Registration 16,663 
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Serial No.: 09/561,272 



.Filing Date: _ 6/7/00 



InvcntorfsV T^1^14 ^^C^V; ^ 



Docket No.:_ 136.147 



Attv.: _JEN;ds 



4/11/02 



The followiiig dated_ 

U,S. Patent Office on Uie DATE STAMPED hereon: 

Amendment 

Issue Fee Transmittal Letter 

Letter 

J-ctier to the OfFicial Draftsperson , 

_Dra«ang5 . Sheets" 



Customer No. 23907 
_ has been received in che 



.Check in the amount of S 



Check No. 



Jnfoimation Disclosure Statement, FonnCs) PTO/SB/08A and PTO/SB/08B & cited referciices 
.rnonty Document 

_OQier_ Rene^^^ p^r^t^ ion u^dfer 37 CFR 1,47 (a) ;DecUratlan & 
— S^f l tPTnpn^ of F^rrq hv Christ laTi Hamon: orleltial cartlfi^d 
^return recelpc postcard ~ 
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Firstar Center, Suite 2000 
777 East Wisconsin Avenue - 
. MilwaBkfce,WI 53202 
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France Telecom ^ Development Division 

France Telecom R&D 



Your correspondent : C.HAMON 

Phone : 02 99 12 48 05 

Ref : FrRDA^AT/VPI/REN/009/Ol/CH 



Mr ALT JALALI 

At Mr HAMID ESMAILFs 
7, JALAL ALLEY 
ASEF STREET 
ZAFARANIEH, VALIASRAV 
TEHERAN (19877) 



Cesson-S6vign6. February 1st, 2001 



Subject : US Patent Application n** 09/581.272 filed December 7. 1998 concerning a "Method of calculating the 
Fast Fourier Transform and the Inverse Fast Fourier Transform", ref CCETT9498; US Declarations; 



Sir, 



Please find enclosed for signature a declaration concerning the US Patent Application you are one of the 
inventors. This signature is indispensable for us, the first and second mails being left unanswered. 

Thank you for sending the executed document back us within three weeks. 

We are looking forward to hearing from you. 



The person in charge for Protection 
in France Telecom R&D Rennes 

C.HAMQN 



France Telecom R&D 

4. rue du Clos-Courtel - BP 59 - 35512 Cesson-S^vign^ Cedex 
Phone : 02 99 12 41 11 Facsimile : 02 99 12 40 98 
Intemational : +33 2 12 41 44 Fax : +33 2 12 40 98 



fr3l1C6 t^i^COi^^ France Tellcom^Dljjj^ 



Votre correspondant : C. Hamon 

Tdi^phone : 02 99 12 48 05 

R^f^rence : FTRDA^ATATI/REN/009/Ol/CH 



Monsieur ALI JALALI 

Chez M. HAMID ESMAILI 
7, JALAL ALLEY 
ASEF STREET 
ZAFARANIEH, VALL\SRAV 
TEHERAN (19877) 



Cesson-Sdvign6, le 1^ fdvrier2001 

OBJET : Demande de brevet aux USA n° 09/581,272 du 07.12. concemant un «proc^d6 de calcul de la 
transform6e de Fourier rapide et de la transform^e de Fourier rapide inverse», ref CCETT9498 ; Pouvoirs aux 
USA 



Monsieur, 

Je vous adresse ci-joint pour signature un pouvoir relatif i la demande de brevet US du brevet dent 
vous ^tes un des inventeurs. Cette signature nous est indispensable, nos premier et deuxiSme envois dtant 
rest^s sans r^ponse. 

Merci de nous renvoyer ce document signe sous trois semaines. 

Dans cette attente, nous vous prions d'agrder. Monsieur, I'expression de nos sentiments distingu^s. 

Le Responsable Protection France Telecom R&D Rennes 




France Telecom R&D 
Site de Rennes 

4, rue du Clos-Courtel - BP 59 - 35512 Cesson-Sevigne Cedex 
Telephone : 02 99 12 41 1 1 - Telecopie : 02 99 12 40 98 
Internationa! : + 33 2 12 41 44 - Fax : +33 2 12 40 98 

SA au capital de 4 098 458 244 EUR - RCS PARIS B 380 129 866 



DECLARATION FOR U^TY OR DESIGN PATENT APP^ATION 

(37 CFR 1.63) 

□ Declaration Submitted with Initial Filing 
OR 

0 Declaration Submitted after Initial Filing (surcharge (37 CFR 1,16(e)) required) 

Attorney Docket Number: 136.147 

First Named Inventor: AH JALALI et al. 

COMPLETE IF KNOWN 

Application Number: 09/581,272 ^" 

Filing Date: June 7, 2000 t 

Group Art Unit: 

Examiner Name: ' " 

As a below named Inventor. I hereby decfare that: 

My residence, post office address, and dtizenship are as stated below next to my name. 

I believe i am the orfglnal, first and sole Inventor (if only one name is listed belowl or an ori«fn^i « . , . 

the specilication of which 
□ is attached hereto 

OR ■ 

0 was filed on June 7. 2000 as United States Application No. 09/581.272 and was amended on June 7. 2000. 
i:r/dtan?amUren'r:^Sr^^^^^^ •^^"'^"'^ °'«^« above-,de„.-ned speci.cation. indlding ihe Cairns, as 

I acknowiedge the duty to disclose infomiation which Is material to the patentability asdefined in 37 CFR 1.56. 
I hereby daim foreign priority benefits under 35 U.S.C. llSfaWdl or aes/hi nf anw ., „ 

certifioate. or 365(a) of any PCT International applicabon whfch de^gS a? I^^^^^^^ «.f„'J'"^.H°"t^ P^**"' 
America. Bsted below and have also identified beick/. by ttx^^ ^tS^r^f^r.^^.^^^,''^^' Elates of 

orofany PCT intema«onai application having a filing da^rbe^Ca1o^°^^^^^^^ 

Prior Foreign Applicatlon(s) p^orftyNot Certfied Copy 

97157-17 ^ Claimed Attached? 

" . Prance December a . 1997 n n v « 

^^"'"''*'> ^^^) (Month/Day/yUr Filed) aVesHNo 



(Number) /Co.,nfrv\ " — tti — TZl^ — r-: ■ ° D Yes □ No 

' (txjuntry) (Month/Day/Year Red) 



□ Add-r Tr^^ '^^^^ (MonthA^ay/YearRled) ° ° ° 

Additional foreign applicaUon numbers are listed on a supplemental prforily data sheet PTO/SB/028 attached hereto: 

hereby Claim the benem under 35 U.S.C. 119(e) of any united States provls^^^^^^ 

(Application Number) {Month/Day/Year Filed) " ° AddiUonal provisional applicaUon 

numbers are listed on a 

7T — - — : supplemental priority data sheet 

(Applicabon Number) (MonUi/OayA'ear Red) PTO/SB/02B attached hereto. 
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DECLARATION m-llty or Design Patent Appllcatl^ 



l^^ti 



Inventors: Ali JALALI, Pierre LeRAY. Dominique LACROIX 
Serial No. 09/581 ,272 



not disclosed in the prior United States o; PCT InSuonal a^SL in ^™„^, 1^^^.°' ^ °' aPPHcatlon Is 
112 I acknowledge the duty to disclose InfoSaSS^S^iSrial to Ca^nS jTd^Rnil.^' °' ^ "-^-C- 

ava..ab.e between the filing date of me prior application and the national or pSSS^li^tE o^tSsl^S^ "'""^ 
U.S. Parent Application or PCT Parent ApplicaBon(s) 4t- - w^uon. 



PCT/FR9aA)2636 December 7. 2000 

(Number) ' (MonthmayA'ear Filed) 



(Patent Number (if appGcable)) 



As a nan^ed Inventor. I hereby appoint the following registered practiBoner(s) to Drosecuta ihi« =,n«i.^H w . ^ 

! ness In the Patent and Trademari< Office connected therewith: Prosecute this application and to transact ali 

NILLES & NILLES, S.C. 
Rrstar Center, Suite 2000 
777 East Wisconsin Avenue 
Milwaukee, Wisconsin 53202-5345 

hereby declare that all statements made herein of my own knowledge are true and th^f »n «f^f^^ * 

)elief are believed to be true; and further that these statements weremLi JS^ Ih^t^ J statements made on InfomiaUon and 
ike so made are punishable by fine or lmpriso^enror"S. unL 18 U^^^ and the 

2opard.2e the validity of the application orany patent issued thereon. statements may 

•ull name of Sole or First Inventor: n a r,o«««« u ^ ^ 

K, '^^''^^^ for this unsigned Inventor 

.iven Name (first & middle nf any]) & Family Name/Surname: Ali JAUVLI 

iventor's Signature: AU JAl^L-X ^ 

Date: ^ 



:esidence (dly, stale, country):. Rennes, France 
Post Office Address: 32 rue Mirabeau 

Bailment K No. 234 
(dty, slate, zip. country): 35700 Rennes. France 



Citizenship: French 
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DECLARATION -Utility or Design Patent Application 



Inventors: Ali JALALI. Pierre LeRAY, Dominique LACROIX 
Serial No. 09/581,272 



Full name of Second Inventor, If any: 

Given Name (first & middle pf anyj) & Family Name/Surname: Pierre LERAY 



□ A peliUon has been filed for this unsigned Inventor 



Inventor's Signature: 
Residence (dty, state, country): Uffre, France 

Post Office Address: LaGrosse Roche 
(dty, state, zip. country): 35340 Uffre, France 




Date: 

Citizenship: French 



Full name of Third Inventor. If any: 



C^' ^n Name (first & middle pf any]) & Family Name/Surname: Dominique LACROIX 



□ A peOtion has been filed for this unsigned inventor 



Dale: lA 1^ / 200n 



Inventor's Signature: ^Vmtnf^uU. /J9CJ^0I ^ 

Residence (dty. stale, country): Rennes, France ^ ^ Citizenship: French 

Post Office Address: 1 0 Cquare Alain Fe i ytt ii t S^jjmfltrt:- L4M\%. Sou^Xi^^r" 




(dty, state, zip, country): OeeoCT Rennes, France 
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ASSIGNMENT ' 

WHEREAS, we, ALI JALALI. PIERB^fcRAY. and DOMINIQUE LACROD^ve made a certain 
new and useful mvention for which we are about to make application for Letters Patent of the United 
States, which application may be identified in the United States Patent Office as Serial 
No. 09/581,272, filed June 7, 2000, and entitled: 

Method of Calculating the Fast Fourier Transform and the Inverse Fast Fourier Transform 

(Executed: , ) 

WHEREAS, France Telecom SA, a company organized and existing under dnd by virtue of the laws of 
France and having its principal place of business at 6, Place d'AJlerayr 75015 Paris. France and 
Telediffusion de France SA, a con^any organized and existing under and by virture of the laws of France 
and having its principal place of business at 10, rue d'Oradour-sur-Glane', 75732 Paris Cedex 15, France 
are.desirous of acquuing the enture interest m and to said invention, said application and the Letters Patent 
to be obtained therefor; 

NOW, THEREFORE, for and in consideration of One DoUar and other good and valuable consideration 
to us in hand paid, the receipt and sufficiency whereof are hereby acknowledged, we have sold assigned 
fnd set over and by these presents do hereby seU, assign, and set over unto the said France Telecom SA 
Telediffusion de France SA, said assignees' legal representatives, successors and assigns, the entire 
nght, title and interest m and to any and all applications including the aforesaid application for Letters 
Patent heretofore and hereafter filed m flie United States or any other country and which may be based in 
whole or in part on said inventions and discoveries, and in and to any and all Letters Patent heretofore or 
hereafter granted by the United States or any other country and which may be based m whole or in part 
on said mventions and discoveries, said assignment including the right to file and prosecute any and all 
such applications and also including the right to sue and recover for any and all infiingements of said 
patents; and request the Commissioner of Patents to issue said Letters Patent to the above-mentioned 
assignee agreeably with the terms of this assignment. 

THE UNDERSIGNED HEREBY GRANT the law firm of NiUes & NUles, S.C. the power to msert in 
this mstmment any further identification which may be necessary or desirable in order to comply with the 
rules of the United States Patent Office for recordation of this document. 

UPON SAID CONSIDERATION, we do hereby covenant for ourselves and our heirs legal 
representatives and assigns that neither we nor any of our said heirs, legal representatives or assig^ have 
or will execute any instrument or perform any act in conflict herewith and that we or our said heirs legal 
•. resentatives and assigns will at aU times do such acts and execute such papers, without expeJise to 
ourselves, as may be necessary or desirable in order to fiilly protect said inventions and discoveries for the 
benefit of said assignee, its successors or assigns, and to otiierwise carry into full force and effect the text 
and interest of this assignment. 



IN WITNESS WHEREOF, we have hereunto set our hands and affixed our seals 
Date: < 



Ali JALALIi 

Date: V g/ Qp/ gg^Qoo 



± (.SEAL) 




Pierre LERAY 



(SEAL) 



Date: ZJRl2Dnn ^T-»^5^^l^ (SEAL) 

Dominique LACRODC '* 



NILLES & NILLES. S.C. 
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FICHE DE DiPdT D'UN ENVOI INTERNATIONAL " ""®'?*^'>* 
^. - AVEC VALEUR DiCUVRiE 

■ A remplir par I'expediteur ^ 



Localite / ^HC^f^ 

Reserve au service 



Pays tf^^ 



signature de I'agent 



RCS Nanhsrre B 354 000 OOO 




□ ^g/'ol standard 
n" de renvoi : 

02755 



. J^^^^'^M^ 



□ colis postal 

ECONOMIQTJE 

valeur declares 
montant : 



□ 



□ 



□ 



l^Ipmpu;; rials on- spcialeV-;;;,Q_;.-; ntnO'V' 



L' env oi m e ntlbrrs^ i C u'^iTriTT frn^ i rtr flfij — • 



□ 



□ 



inscrlt sur CCP 



Date et signature- 



X Code ipbist^ et;- Ipcalite p-;-; 






Telecom 



Your correspondent : C.HAMON 

Phone : 02 99 12 48 05 

Ref : FTRD/DSTATI/REN/l 13/00/CH 



Development Division 
France Telecom R&D 



MrALI JALALI 

At Mr HAMID ESMAILI's 
7, JALAL ALLEY 
ASEF STREET 
ZAFARANIEH, VALIASRAV 
TEHERAN (19877) 



Cesson-S^vign^. September 28. 2000 



Subj^t : US Patent AppUcation n° 09/581,272 filed December 7. 1998 concerning a "Method of calculating the 
Fast Founer Transform and the Inverse Fast Fourier Transform", ref CCETT9498; US Declarations; 



Sir, 



Please find enclosed for signature a declaration concerning the US Patent Application you are one of the 
mventors. This signature is indispensable for us. the first mail being left unanswered. 

Thank you for sending the executed document back us within three weeks. 

We are looking forward to hearing fi^om you. 



The person in charge for Protection 
in France Telecom R&D Rennes 

C.HAMON 



France Telecom R&D 

4. rue du Clos-Courtel - BP 59 - 35512 Cesson-S^vign6 Cedex 
Phone : 02 99 12 41 1 1 Facsimile : 02 99 12 40 98 
International : +33 2 12 41 44 Fax : +33 2 12 40 98 



t 



f ranee teleCO'^ BrancheO^veloppem 



France T6l§com R&D 



Votre correspondant : C. Hamon 

T6I6phone : 02 99 12 48 05 

R^f^rence : FTRD/DSTVATI/REN/l 13/00/CH 



Monsieur ALI JALALI 

Chez M. HAMID ESMAILI 
7, JALAL ALLEY 
ASEF STREET 
ZAFARANIEH, VALIASRAV 
TEHERAN (19877) 



Cesson-S6vign6. ie 28 septembre 2000 

OB JET : Demande de brevet aux USA ii° 09/581,272 du 07.12. concemant un «proc^d6 de calcul de la 
transform^e de Fourier rapide et de la trasnform^e de Fourier rapide inverse», ref CCETT9498 ; Pouvoirs aux 
LISA 



Monsieur, 

Je vous adresse ci-joint pour signature un pouvoir relatif A la demande de brevet US du brevet dont 
vous etes un des inventeurs. Cette signature nous est indispensable, notre premier envoi 6tant restd sans 
rdponse. 

Merci de nous renvoyer ce document signd sous trois semaines. 

Dans cette attente, nous vous prions d'agr^er. Monsieur, I'expression de nos sentiments distingu^s. 



Le Responsable Protection France T616com R&D Rennes 




Christian HAMON 



France Telecom R&D 
Site de Rennes 

4, rue du Clos-Courtel - BP 59 - 35512 Cesson-Sevign6 Cedex 
Telephone r 02 99 12 41 1 1 - T§!6copie : 02 99 12 40 98 
International : + 33 2 12 41 44 - Fax : +3r3 2 12 40 98 

SA au capital de 4 098 458 244 EUR - RCS PARIS B 380 129 866 



WemKA I lUN hbWWIlTY OR DESIGN PATENT APplB 

(37CFR1.63) 



fATION 



□ Declaration Submitted with Initial Filing 
OR 

0 Declaration Submitted after Initial Filing (surcharge (37 CFR 1 .1 6(e)) required) 

Attorney Docket Number: 136.147 

First Named Inventor: Ali JALALI et al. 

COMPLBTE IF KNOWN 

Application Number: 09/581,272 

Filing Date: June 7. 2000 ^ i. 

Group Art Unit: 

Examiner Name: . \ \ 

As a below named Inventor, I hereby declare that: 

My residence, post office address, and citizenship are as stated below next to my name, 

I believe I am the original, first and sole Inventor flf only one name is listed below) or an original, first and joint inventor fif olurai 
names are listed belc^) of ttie subject matter which Is claimed and for which a patent is sought on t^iernvSoren?^^^^^^ ^ ' 

Method ofCatculatmg the Fast Fourier Transform and the inverse Fast Fourier Transfomi 
the specification of which 
□ is attached hereto 
OR 

0 was filed on June 7. 2000 as United States Application No. 09/581 .272 and was amended on June 7. 2000. 

^"^^"^ above-idenUfied speancaUon. including the Cairns, as 

I acknowledge the duty to disclose infonnation which is material to the patentability as defined in 37 CFR 1 .56. 

'c^Ste"'TS7Z^^^ ''^^^^1} appiicaUon(s) for patent or inventor's 

A^n^^Z. , • °^ intemabonai application which designated at least one country other than the United States of 

^of^J p''?.^''.^'^^^"^ .^^^^ '^^"""^ ^^'"^^ the box. any foreign application for pLr^^^^^^ 

or of any PCX international application having a filing date before that of the application on which priority Is d^^^^ certificate. 

Prior Foreign Application(s) 

^^"^^Z France Decembers, 1997 

(Number) (Country) (Month/Day/Year Filed) 



(Number) (Country) (Month^ayA^ear Filed) 



Priority Not 


Certified Copy 


Claimed 


Attached? 


□ 


□ Yes 0 No 


□ 


□ Yes n No 


□ 


□ Yes □ No 



(Number) (Country) (Month/Day/Year Filed) 

□ Additional foreign application numbers are listed on a supplemental priority data sheet PTO/SB/02B attached hereto: 
I hereby claim the benefit under 35 U.S.C. 119(e) of any United States provisional application(s) listed below. 

(Application Number) (MontiVDavA'ear Filed) ° Additional provisional application 

' numbers are listed on a 

supplemental priority data sheet 



(Application Number) (Montii/DayrTear Filed) PTO/SB/02B attached hereto. 
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DECLARATIOI^^Jtility or Design Patent Applicatil^ 

Inventors: Ali JALALI. Pierre LeRAY. Dominique LACROIX 
Serial No. 09/581,272 

I hereby claim the benem under 35 U.S.C. 120 of any United Stales appllcation(s), or 365(c) of any PCT international application 
designating the United. States of America, listed below and. Insofar as the subject matter of each of the dalms of this aDoncaflon is 
not disclosed in the prior United States or PCT IntemaUonal application In the manner provided by the first paragraph of 35 U S C 
1 12, 1 acknowledge the duty to disclose information which is material to the patentability as defined In 37 CFR 1 56 which became 
available between the filing date of the prior application and the national or PCT international fillna date of this applicatipn. 
U.S. Parent Application or PCT Parent Application(s) 

PCT/FR98/02636 December 7. 2000 , ^ 

(Number) (MonthA>ay/Year Filed) (Patent Number 0f applicable)) 



(Number) (MonthAJay/Year Red) (Patent Number Of applicable)) 

O Additional U.S. or PCT International application numbers are listed on a supplemental prfority data sheet PTO/SB/02B attached 

ri6r6i0> 

As a named Inventor, I hereby appoint the fonowing registered practiUoner(s) to prosecute this application and to transact all 
business in the Patent and Trademark Office connected therewith: uansaci an 

James E. Nilles. Reg. No. 16,663 Matthew C. Loppnow, Reg. No 45 314 

Andrew J. Nilles. Reg. No. 31.786 Stephen Michael Patton, Reg. No. 36 235 

Jay G. Durst, Reg. No. 41,723 Jerome D. Drabiak, Reg. No. 31 Oil ' 

Usa M. Gehrke. Reg, No. 38,888 Lisa A. Brzyckf, Reg. No. 40,926 
Thaddeus C. Stankowski. Reg. No. 45.522 

Direct all telephone calls to James E. Nilles at telephone number (414) 276-0977. facsimile number (414) 276-0982 
Direct all coaespondence to: James E, Nilles ^^o^. 

NILLES & NILLES, S.C. 

Firstar Center, Suite 2000 

777 East Wisconsin Avenue 

Milwaukee. Wisconsin 53202-5345 

I hereby dedare that all statements made herein of my own knowledge are true and that all statements made on information and 
belief are beheved to be tme; and further that these statements were made v^th the knowledge that v^llM fel^ st^terS 

Lop"^^^^^^^^ '?"f, or both, under 18 U.S.C. 1001 and Lt such v^ilfulfelse S^^^ 

jeoparaize the validity of the application or any patent Issued thereon. 

Full name of Sole or First Inventor: □ a peUtion has been filed for this unsigned inventor 

Given Name (first & middle [if any]) & Family Name/Surname: AJi JALALI 



Inventor's Signature: ^ Li 'ZP^U\L-\ ^ 



Date: ^ 



Residence (city, state, country): Rennes, France Citizenship: French 

Post Office Address: 32 rue Mirabeau 

Batiment K No. 234 
(city, state, zip. country): 35700 Rennes, France 
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DECLARATIO^^Jtility or Design Patent ApplicatJ 



Inventors: Ali JALALI. Pien-e LeRAY. Dominique LACROIX 
Serial No. 09/581,272 



Full name of Second Inventor, if any; 
Given Name (first & middle pf any]) & Family Name/Surname: Pierre LERAY 



□ A petition has been filed for this unsigned inventor 



Inventor's Signature: ^^ (^rre 
Residence (city, state, country): Uffre, France 

Post Office Address: LaGrosse Roche 
(city, state, zip. country): 35340 Uffre. France 




Date: 

Citizenship; French 



Full name of Third Inventor, if any: 



Given Name (first & middle (if any]) & Family Name/Surname: Dominique LACROIX 



□ A peUtion has been filed for this unsigned Inventor 



Date: lAI^I^OOn 



Inventor's Signature: ^VmfntaXAJL /j^CJ^Ol ^ 

Residence (city, state, country): Rennes. France ^ ^ Citizenship: French 

Post Office Address: 1 8 Squa r e Alain re i ye ii l Sa S<|jwOLrt=. Z^fl ^oSxx^\£^ 




(city, state, zip. country): 36600 Rennes. France 
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ASSIGNMENT 

WHEREAS we. ALI JALALI. PIE^fc^ERAY. and DOMINIQUE LACRolHave made I certain 
new and useful invention for which we are about to malce application for Letters Patent of the United 
States, which application may be identified in the United States Patent Office as Serial 
No. 09/581,272. filed June 7, 2000, and entitled: 

Method of Calculating the Fast Fourier Transform and the Inverse Fast Fourier Transform 

(Executed: ) 

WHEREAS, France Telecom SA, a company organized and existing under tod by virtue of the laws of 
France and having its principal place of business at 6. Place d'Allerayr 75015 Paris. France and 
Teledififusion de France SA, a company organized and existing under and by virture of the laws of France 
and having its principal place of business at 10, rue d'Oradour-sur-Glane', 75732 Paris Cedex 15 France 
are desirous of acquiring the entire interest in and to said invention, said application and the Letters Patent 
to be obtained therefor; 

NOW, THEREFORE, for and in consideration of One DoUar and other good and valuable consideration 
to us m hand paid, the receipt and sufficiency whereof are hereby acknowledged, we have sold assigned 
and set over and by these presents do hereby sell, assign, and set over unto the said France Telecom SA 
and Telediffiision de France SA. said assignees' legal representatives, successors and assigns, the entire 
nght. title and mterest m and to any and all applications mcluding the aforesaid application for Letters 
Patent heretofore and hereafter filed in the United States or any other country and which may be based m 
whole or in part on said inventions and discoveries, and in and to any and aU Letters Patent heretofore or 
hereafter granted by die United States or any other country and which may be based m whole or in part 
on said inventions and discoveries, said assignment including the right to file and prosecute any and all 
such applications and also including the right to sue and recover for any and all infringements of said 
patents; and request die Commissioner of Patents to issue said Letters Patent to the above-mentioned 
assignee agreeably with die terms of this assignment. 

THE UNDERSIGNED HEREBY GRANT the law firm of NUles & NiUes, S.C. the power to insert in 
this mstrument any further idendfication which may be necessary or desirable in order to comply with die 
rules of die United States Patent Office for recordation of diis document. 

UPON SAID CONSIDERATION, we do hereby covenant for ourselves and our heirs legal 
representatives and assigns that neither we nor any of our said heirs, legal representatives or assigiis have 
or will execute any instrument or perform any act in conflict herewifli and that we or our said heirs legal 
representatives and assigns will at all times do such acts and execute such papers, widiout expeJise to 
ourselves as may be necessary or desirable in order to fully protect said inventions and discoveries for die 
oenetit of said assignee, its successors or assigns, and to otherwise carry into full force and effect die text 
and mterest of this assignment. 



IN WITNESS WHEREOF, we have hereunto set our hands and affixed our seals. 
Date: < 



Date: ^^/oS'/^ooo 

Pierre LERAY 




(SEAL) 
(SEAL) 



Dominique LACROK 

NILLES & NILLES S C 

O.^. , \\NTSERVER\COMPANY\Oau\CUEmM3tfV147\ASSIGNMENT.doc 



UAbSINt:! BALLUT-SCHMIT 

Brevets - Marques - Dessins - Modeles - Contrats - Consultations - Litiges 
4, rue Generqjgli phe - BoTte Posta leSSIS - 56^^ 

" '■ ^^^ ^' ^^''^^^ - ll^cop!e02^76^ 

http://www.balIbt-schmitcom - e.mall:info@baHot-schmit.com 




FRANCE TELECOM - CNET 

DSTVAO^I - Service Gestion 

38-40, rue du General Leclerc 

92794 ISSYMOULINEAUX CEDEX 9 

A Vattention deM. Didier LEMOYNE 



Copie 

Vos ref 
Dossier 
Lettre 



Lorient, le 10 aout 2000 

Mr HAMON (site Rennes). 

9498/CNET2966 
013770 US/pc 
BH/ml/00-0760. 



Demande de brevet aux Etats-Unis n"" 09/581 272 



Depos6e 
Titulaire 
Inventeur 
Priorite 
Titre 



07/12/1998 

FRANCE TELECOM (SA) et TELEDIFFUSION DE FRANCE (SA) 
JALALI Ali - LERAY Pierre - LACROIX Dominique 
France - n° 97 15737 - 08/12/1997 

PROCEDE DE CALCUL DE LA TRANSFORMEE DE FOURIER 
RAPIDE ET DE LA TRANSFORMEE DE FOURIER RAPIDE 
INVERSE. 



Cher Monsieur, 

Nous nous permettons de vous rappeler notre courrier du 05/07/2000 
concemant le passage en phase nationale aux Etats-Unis de la demande 
mternationale PCT/FR98/02636. 

Comme nous vous le disions alors, afm de regulariser cette affaire, il est 
indispensable de deposer un pouvoir et une cession signes par les inventeurs a 
savoir Madame Dominique LACROIX et Messieurs Ali JALALI et Pierre 
LERAY. 

Ces documents ne nous etant toujours pas parvenus, le Patent Office americain 
vient de nous transmettre une notification fixant le delai de presentation desdits 
documents au 26 aout prochain. 

Nous vous demandons done de nous les retoumer le plus rapidement possible 
afin qu'ils puissent etre deposes au Patent Office dans les delais impartis. 

Dans Tattente de vous lire, nous vous prions d*agreer. Monsieur, I'expression 
de nos sentiments les meilleurs. 



^ Bertrand HAYS 



SOCIETE DE CONSEILS EN PROPRIETE INDUSTRIELLE 
S.A. au capital de 250.000 F - 7. rue Le Sueur - 751 16 PARIS 
RCS Paris B 341 391 308 - RCS Lorient 90 B 566 - APE 741 A 

CONDITIONS DE RfeGLEMENT Suivant tusago des professions libirales. les ordres ne sent exdcutds qu'apfds le rdglement de leur moniant. 



Avis & coller p»x le guic&eil^r 




S9I SBd SO^^iJK 



France Telecom 




Development Division 
France Telecom R&D 



Your correspondent : C. Hamon 



Phone: 02 99 12 48 05 

Ref : FTRDA^AT/PI/REN/062/Ol/CH 



MrALI JALALI 



At Mr HAMID ESMA-ILrs 
7, JALAL ALLEY 
ASEF STREET 

ZAFARANIEH, VALIASRAV 
TEHERAN (19877) 



Cesson-Sevigne, June 26, 2001 



Subject : US Patent Application n** 09/581,272 filed on December 7, 1998 conceming a "Method of calculating 
the Fast Fourier Transform and the Inverse Fast Fourier Transform", ref CCETT9498; US Declarations. 



Please find enclosed for signature a declaration conceming the US Patent Application you are one of the 
inventors. This signature is indispensable for us, the first and second and third mails being left unanswered. 

We enclose a copy of the US patent application as well as the Preliminary Amendment. 

Thank you for sending the executed document back us within three weeks. 

We are looking forward to hearing from you. 



France Telecom R&D 
Office of Rennes 

4, rue du Clos-Courtel - BP 59 - 35512 Cesson-Sevigne Cedex 
Phone : 02 99 12 41 11 Facsimile : 02 99 12 40 98 
International : +33 2 12 41 44 Fax : +33 2 12 40 98 



Sir, 



The person in charge for Protection 
in France Telecom R&D Rennes 



Christian HAMON 




f ranee teleco 



Branche DSveloppeme 
France Telecom R&D 



Votre correspondant : C. Hamon 

T^Idphone : 02 99 12 48 05 

R^fdrence : FTRDA^AT/PI/R£N/062/01/CH 



Monsieur ALI JALALI 

Chez M. HAMID ESMAILI 
7, JALAL ALLEY 
ASEF STREET 
ZAFARANIEH, VALL\SRAV 
TEHERAN (19877) 



Cesson-S^vignd, le26 juin 2001 



^5oL.e dfpovSetrJid^T^^^^ 07^12. concen^ant ua «proc.d. de calcul de la 

USA transform^e de Founer rapide mverse». ref CCETT9498 ; Pouypirs aux 



Monsieur, 

Je vous joins une copie de la demande de brevet aux USA ainsl que ramendementpr^Iiminaire . 
Merc! de nous renvoyer ce document sign^ sous trois semaines. 

Dans cette attente, nous vous prions d'agr^er. Monsieur, I'expression de nos sentiments distinguds. 

Le Responsable Protection France T^l^com R&D Rennes 




Christian HAMON 



France Telecpra.R&D.- 
Site de Rennes 



4 ue du Clos-Courtel - BP 59 - 35512 Cesson-Sevigne Cedex 
Telephone : 02 99 12 41 11 - Telecopie : 02 99 12 40 98 
International : + 33 2 12 41 44 - Fax : +33 2 12 40 98 
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LJ XV IS Da RBCTPTIOi* i^ UW ENVOI INTERKATIONAL . 

□ AVIS DE PAIEMEirr / INSCRIPTION DUN MANDAT SUB UN CCP 



Bureau de ddpdt 



Date de dSpdt 



A rempUr par I'Q t^M^r tS'T^E ST' , ^ 



Service courrler 



n envoi standard 
□ PRIORITAIRE 
S n'de renvoi 



postal LI Uvre LJ sac M 



CW 07 
Service des Postes 

PRIORnAIRE / PAR AVION 



r 



Timbre du bunau 
renvoyant I'avla 



□ coUs 

□ fiCONOMIQDE 

□ vaJeur ddclar6e 
montanb : 



•^H'^^4<envoyer & 



Services flnahclers 



□ mandatn" 



montant : 



L. 



_l 



A completer & destination 



L'envol mentionnS ci-dessus a 6t6 dOment : 

Q remls Ell pay^ Q Inscrit sur CCP 



Date 6t signature * 



Norn ou raison soclale 



Rue et n** 



Code postal et locality 



FRANCE 



' * Cet avis pourra etre sl^d par le desUnataire ou, si les r^^emeats du pays de destination 
j'ls pMvoient, par une autre personne autorlsSe ou par I'a^ent du bureau de desUnatioa. 



608S33 B • C]r07 - VI - 9923147470012300 - D17 • B08833 



JLPOSTK 



estinataire 



FICME DE DiPOT D'UN ENVOI INTERNATIONAL 



- RECOMMANDi 

- AVEC VALEUR D^CLARIE 



—————— A remplir par I ' expe d i t eu r 



Localite 



Pays T/Ph/ 



Reserve au service ■ 




date N° 


1 Prix 1 vd 1 


Nature 


signature de I'agent 












o 














Poids 

1 1 
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° °°°'oR "°" '""tal Filing 

0 Dec.a.a»o„ Su,m,,,ed a«er l„l„a, Filing (s„.cha.oe (37.CPR ,.16(e„ required, 

Attorney Docket Number: 136.147 

First Named Inventor; Ali JALALI et af 

COMPLETE IF Kf/OVW J/*»-AU et al. . , 

.Appfication Number: 09/581 272 O- 

Filing Date: June 7.' 2000 • 

Group Art Unit: , >■ 

Examiner Name: '■ 

As a below named Inventor. I hereby declare lhaf 

My residence, post office address, and cffeenship are as slated below next t 

« le specification of which 'nverse f^ast Founer Transfom 

O Is attached hereto 
OR 

0 ^^«'^<^°"J"ne7,2000asUniWStatesAppncalIonNo 09/58^^ 

(Month/Day/YearFBed) ' OYesHNo 



(Number) t:; r-r— m 

^^""'^J (Month/Day/YearRled) DYesDNo 

(Number) ^ — p. 

Additional foreign application numbers arSed'^l . , (Month/Uay/Year R.ed) " ° ° ° N" 

^ '""P'^'"®"'^' P'tority data Sheet PTO/Sa/02B=.H . 
erebydaimthebenefitunder35USC , ™SB/028 attached hereto: 

/--"^"•^•^• ^^9(e)of any United States pro.-s,ona,.applica.ion(s),istedbe,ow. 

(Application Number) ' — n 

(Month/Day/Year Filed) ° '^^'"^^^"a' Pravfsional applicabon 

numbers are fisted on a 



(Application Number) 7 M^^m/r. ^^iPP'^"^^"^' Priority data sheet 

(Monin/uayA-ear Fifed) " PTO/SB/02B attached hereto 
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W^esign Paten. Appi^oatoP"'' ^ 

I hereby dafm Uie benefit under ssutf- ^-tn , 

lesIgnaSng the United Sta^S^±V, ° h """^'^ ^^^^^ appflcat;on(s). or ass^ci nf o^. 

f12 I acknowledge the duty to disctosahSJon^Si f ^PP"«''*«" 'tanner prowJSbfth AS ™ °^ "^'^ appBcatfon b 
-vailable between the filing date of Uie pri2?SS^L?A'^^^^^ "'^ patenlabffij as dSnedt i^^co^^"^'' ^-S-C 

'-ParentAppKcat^norPCTPareneC^^^^^^^^^ 
PCT/FRQn/ noet-ic 

(Month/uay/Year Filed) " .p,. ^ 

. , ■■ (Patent Number (.fappficablej) ~ 

(Number) pr; — jrrr: ■ 

James E. Nilles. Reg. No. 16.663 ^ . 

^drew J.Nilles. Reg. NO. 31.786- S^m," No. 45.314 

Jay G Durst. Reg..No. 41.723 - jjf^^"r!^'^^^'^^««^- ReS-.N^^ 
Usa M.Gehri<e. Reg. No. 38.888 SfT p^' No. 31.011 

Thaddeus C. StankowsW. Reg. No. 45.522 Brzyckl. Reg. No. 40,926 

NfLLES& NILLES. S.C. ^276^98^ 
Rrstar Center, Suite 2000 
777 East Wisconsin Avenue 
Milwaukee, Wisconsin 53202-5345 

name of Sole or First Inventor ^ 
Name (first & middle pf any], & Family Name/Surname: Ali MUVLI " ^ ""^ "'"'S-^ ""-^n'o^ 

M's Signature: J4 tyj- £_i X 

'"ce (dty. state, country):. Rennes. France ^ 



Post Office Address: 32 rue Mirabeau Citizenship: French 

Bailment K No. 234 
(City, state, rtp. country): 35700 Rennes. France 
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DECLARATION -^ty or Design Patent Appllcatio 



Inventors: Ali JALALI, Pierre LeRAY, Dominique LACROIX 
Serial No, 09/581,272 



Fuir name of Second Inventor, If any: 

GiVeri Name (first & middle [if any]) & Family Name/Surname: Pien-e LERAY 



:nv6ntcf's Signature: _ 
Residence (city, state, country): Uffre, France 

Post Office Address: LaGrosse Roche 
(dty, state/zip, country): 35340 Uffre. France 




O A peb-tion has been fiied for this unsigned Inventor 



^ Date: -Y 6/6 S /<gg?Oo 

Citizenship: French 



•ull name of Third Inventor. If any: 



Siven Name (first & middle (if any]) & Family Name/Surname: Dominique U^^CROIX 



n ApeUtion has been filed for this unsigned inventor 



iven 



tor's Signature: 3^X)mfr\t^bUJL /^CJ^Ol Jk' 




esidence (dty, state, country): Rennes, France 

Post Office Address: 1 0 Cquar i i Aldiii reiye i it Sa SojuiSLn^ 6ov^txrtfj€r- 

V. state. 7Tn miinfn/\* '%CA/WT Dannie. c»-^»^» I Q 



Date:. lA/^I^Onn 
Citizenship: French 



(city, state, zip, country): 36eoa Rennes, France 
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ASSIGNMENT ^ Zt) 

Me,^ Of Calcumni the Fas, Fourier Transform an, rte Inverse Fas, FouHer Transfonn 



(Executed; 



are desirous of acquiring the entire interest in and to said invenh^n cIm V • ^^^^ ^^"^6. 
to be obtained therefor; invention, said apphcation and the Letters Patent 

NOW. THEREFORE, for and in consideration of One DoUar anH nfh. ^ . 
CO us in hand paid, the receipt and sufficiency wSreoTaS £re^ f ^ f ^ ^"'""^^^ consideration 

and set over and by these presents do hereby ITa Z have sold, assigned. 

r 1 Teledifiusion de Fran^ SA. said ^Z^\^?rco^in^tT' ""'^ ^^^^ telecom SA 

r^ht. title and interest in and to any ^ aHopSalTi^SSr^^ ^^s^rs assigns, the entire 
Patent heretofore snd hereafter filed in the l^L^ sStes^r^v o^^^^ application for Letters 

«vhole or mpatt on^aid inventions and disced W S ^i^l^^ f^t"^ be liased in 
lereafter granted by the United States or any oS Sun^^l ^\ ^ heretofore or 

>n said inventions Ld discoveries. siS^Sr^entSri S " °^ - 

«ch apphcations and also mcluding the right to sue and Jernv^^ .^"^ prosecute any and all 

.atenc: and request the Conunissioner o^pIt^L to ^le^aif^^^ infringements of said 

ssignee agreeably with the terras of this assignment. " *° above-mentioned 

HE UNDERSIGNED HEREBY GRANT the law firm of Nilles & Nilles ^ r 

L-frSsrr=sSfj-^^^^ 

^;in^rs^?S legal 
■ will execute any mstrument or perform any a7in coSSlrewil^^^^^^^ ^^^f^^^^^^^ °^ ^«gns have 
presentatives and assigns will at all times do such aS Id ex^^^ °' 

Ives, as may be necessary or desirable in ord^ to tuv^ro^TTri • ^'^"^ to 

nefit of said assignee, its successors or assTg^ and to oAr.SS a'^d discoveries for the 

d interest of this assignment. ^ ' otherwise carry mto full force and effect the text 

WITNESS WHEREOF, we have hereunto set our hands and affixed our seals, 
te: < 

^ (SEAL) 




Fierre LERAY (SEAL) 

'■■-Sdnj2QQQ_ ^^^JS^^^L^ (SEAL, 

Dommique LACROIX ■ Ioeal,; 

'LES & NILLES. S.C. 



^^NTSERVER^COMPANY^D«.^CL.EmM3^^M7VASS.CNMENT., 
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- IN THE UNITED ST^S PATENT AND TRADEMARK OFFICE 
. National Phase of PCT/FR98/02636 

International Filing Date: 07 December 1998 

Inventors: Ali JALALI. Pierre LERAY. and Dominique LACRODC 

Priority: French Application No. 97 15737 
Filed 08 December 1997 

Attorney Docket 136.147 

PRELIMINARY AMENDMENT 

DO/EO/US 

Assistant Commissioner for Patents 
Washington, D.C. 20231 

Sir: 

This Prelhnmary Amendment is directed to a new U.S. application as identified above. 
Please enter this preliminary amendment prior to calculating the fees. 

Please substitute the attached claims 1-42. which incorporate amendments made to 
claim2 during international prehmmary examination, for pages 36-51 containing claims 1-42. 
Please use the substitute clauns for examination purposes. 

Please amend the specification, substitute claims, and Abstract as follows: 

IN THE SPECIFICATION 

Page 1, after the title insert the heading - BACKGROUND OF THE 
INVENTION -; and the subheading - 1. Field of the Invention 
between lines 5 and 6. insert the subheading 

- 2. Description of the Related Art ~; 



liminary Amendment - Ali JALALI et 

thod of Calculating the Fast Fourier Transform... 

;e2 

Page 7. between line. 28 and 29, insert the heading - OBJECTS AND SUMMARY 
'THE INVENTION-; 

■ Page 10, between lines 26 and 27, inset, the heading - BRIEF DESCRIPTION OF 
SE DRAWINGS -; 

Page 13. between lines-S and 9. insert the heading - DESCRffTION OF THE 

^ 

REFERRED EMBODIMENTS --. 

IN THE CLAIMS. As Amen ded Under Article 34 
Claim 3 , lines 1 and 2, cancel "or 2" ; 

Claim 5. lines 2 and 3. cancel "in mrn dependent on claim 3. in turn dependent on 
:laiml."; 

Claim 6. lines 1. 2 and 3. cancel "in turn dependent on claim 3. in turn dependent on 
;laim 1, or according to claim 5," 

Claim 23. lines 1 and 2, cancel "or 20"; 
Claim 25, line 2, cancel "or 24"; 
Claim 26. line 2, cancel "or 25"; 

Claim 29. lines 2 and 3. cancel "in turn dependent on claim 3. in mrn dependent on 

claim 2," 

Claim 31, line 2, cancel "or 30" ; 
Claim 35, line 2, cancel "or 34"; 
Claim 39, line 2. cancel "or 38"; 
Claim 40, line 2, cancel "or 39". 



Pfeliminary Amendment - Ali JALALI 
Method of Calculating the Fast Fourier 
Pages 



IN THE ABSTRACT 
Please substitute the attached Abstract of the Disclosure for the Abstract as filed. 

REMARKS 

This application has been amended to incorporate the modifications made to the 
international application under Article 34 which include changes made to claim 2. The 
application is further amended to insert headings in the specification, elLcninate the multiple 
clarni dependencies, and to conform the Abstract in accordance with U.S. Patent Office 
practice. Entry of the amendnients and early consideration and allowance are respectfully 
requested. 

Respectfully submitted, 

^ James E. Nilles 

Registration No. 16,663 

Dated: June 7, 2000 

NILLES «& NILLES, S.C. 
"^7 East Wisconsin Avenue, Suite 2000 
Milwaukee, WI 53202 
Telephone (414) 276-0977 
FacsimUe (414) 276-0982 

ds G:\Data\CLIENT\136\147\PreLAmend.doc 
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METHOD OF CALCULATING THE FAST FOURIER TRANSFORM 
AND THE INVERSE FAST FOURIER TRANSFORM 

Abstract of the Disclosure 

A method of calculating the fast Fourier transform or the inverse fast Fourier transform 

of a series of N real samples x(i^, with N power of two, operatmg according to a time 

interleaving algorithm and providing the sample series X(n) in ascending order to index n and 

usmg limited calculating and storage means. A method of calculating the fast Fourier 

transform or the inverse fast Fourier transform of a series of N conjugated complex samples 

X(n), with N power of two, operating according to a frequency interleaving algorithm. 



• 3S • 



AMENDED CLAIMS 

1. A method of calculating the fast Fourier 
transforn or the inverse fast Fourier transform of a 
digital signal defined by a series of N real starting 
samples x (n) , with N a power of two and n e [0..N-1], 
comprising successive transformation steps (2) for 
transforming ioput samples into output samples, all the 
transformation steps being performed by means of a 
single set of butterflies with several inputs and 
several outputs, the operating mode of which is 
modified selectively in each transformation step, the 
input and output samples of each transformation step 
being stored in a storage memory, a series of N output 
samples y(n) representative of the fast Fourier 
transform or the inverse fast Fourier transform of the 
output samples x (n) being provided in the last 
transformation step, 

characterized in that output samples y(n) are 

real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x(n) processed in the first transformation step 
are classified in bit-reversed order of their index n, 
output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these oytput samples being defined by the following 
relations : 

y(0) = Re[X(0) ] 



y(n) = Re[X((n+l)/2)] for n being odd and 

different from N-1 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = RetX (N/2) ] 
where samples X(n), with n e [0...N-1], designate the 
complex samples of the series corresponding to the fast 
or inverse f^st Fourier transform of the starting 
sample series x(n). 

2. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N complex samples 
X(n) conjugated by pairs, characterized in that the 
calculation is done on a series of N real starting 
samples y(n) representative of the series of complex 
samples X(n), with N power of two and n e [0..N-1], 
the starting samples y(n) being defined as follows: 

y (0) = Re[X{0) ] 

y(n) = Re[X((n+l)/2)] for n being odd and 

different from N-1 

y(n) = Im(X{n/2)] for n being even and 

different from 0 

y(N-l) = Re [X (N/2) ] 
in that this method comprises successive transformation 
steps for transforming input samples into output 
samples, a series of N real output samples x(n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 



is modified selectively in each transformation step, 
the input and output samples of each transformation 
step being stored in a storage memory, 

and in that the output samples of a butterfly 
5 replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x(n) are provided in the last 

10 transformation step in bit-reversed order of index n. 

3- The calculation method according to claim 1 or 
2, characterized in that, in each transformation step, 
each butterfly transforms input sample pairs, the ranks 
of the input samples of the same pair within the series 

15 of input samples of said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 
butterfly. 

4. The calculation method according to claim 3, 
characterized in that it comprises p-1 transformation 
steps Ep with vi=log2(N) and p e [0..U-2]. 



20 



5. The calculation method according to claim 4, in 
turn dependent on claim 3, in turn dependent on claim 
25 1, characterized in further comprising: 

- a preliminary step of modifying the sequence of 
the starting samples x(n) ranked in ascending order of 
index n and showing them in bit-reversed order of index 
n in the first transformation step, and 

- a final step of processing the series of output 
samples y(n) and providing a series of N complex 
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conjugated samples X(n) corresponding to the fast or 
the inverse fast Fourier transform of the series of 
starting samples x(n), 

6, The calculation method of claim 4, in turn 
5 dependent on claim 3, in turn dependent on claim 1, or 
according to claim 5, characterized in that, in each 
transformation step Ep, butterflies are distributed 
among N/2P"^^ calculation blocks / 

in that each calculation block has a peripheral 

10 butterfly and/or 2^-1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank Z^^^o,, 2P^^a+2P'-^-l, T^^^cl^Z^^^ , 
. 2'^^^a+2'^'^^-l into output samples of the same rank, 

15 and in that the internal rank t butterfly of the 

rank a calculation block in transformation step Ep 
transforms the input samples of rank 2'^*^a+2T+l, 
2P^2a+2t+2, 2P''2a+2'^*2-2x-3, 2P^^a+2*^^2-2x-2 into output 
samples of the same rank, with p>l . 

20 7. The calculation method according to claim 6, 

characterized in that each butterfly is assigned a 
coefficient W^, whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j(2ns/N) ^ith s G [0.. N/4-1] for a fast Fourier transform 

25 and is equal to e^*^"^^^* with s e [0. •N/4-1] for an 
inverse fast Fourier transform. 

8. Calculation method according to claim 7, 
characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 

30 is assigned coefficient with 5 = (x+1) . (N/2P*^) . 
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9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - four inputs for receiving input samples and 

four outputs for providing output samples, 

- four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

10 in order to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 

15 additional inputs. 

10. The calculation method according to claim 9, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

20 in that the permutation signal is 0 for the even 

rank calculation blocks, including rank 0, and 1 for 
the other ones. 

11. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 

25 calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies. 

12. The calculation method according to claim 11, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 

30 transformation step, and otherwise 0. 
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13. The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A4-j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: si = el + e2 

s2 = el " e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 
• s4 = (el + e2) - (e3 + e4) 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si == el + A.e3 - B,e4 
s2 = e2 + B.e3 + A,e4 
s3 = el - A.e3 + B»e4 
s4 = -e2 + B.e3 + A.e4 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 



si = 



el - 



A,e3 + B.e4 



s2 = - 



e2 + B.e3 + A,e4 



s3 = 



el - 



A.e3 " 



B.e4 



s4 = e2 + B.e3 + A.e4 
14." The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 
calculation block comprises: 



2^-1 internal butterflies and a peripheral 
butterfly for the even values of index p as well as for 
the last transformation step if \x is even, and 
2^-1 internal butterflies, otherwise, 
5 15. The calculation method according to claim 13, 

characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 
transformation step with ]i being odd, and otherwise 0. 

16. The calculation method according to claim 15, 
10 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if primary mode, secondary mode and permutation 
15 signals are 0; 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 
20 2) if the primary mode signal is 0 and the 

secondary mode signal is 1: 
si = el + e4 
s2 = e2 
s3 = e3 
25 s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
's2 = el - e2 
30 s3 = e4 - e3 

s4 = el + e2 + e3 + e4 




4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A,e4 
s3 = el - A.eS + B.e4 
s4 = -e2 + B.eS + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1 : 

si = e^^l - A.eS + B.e4 
s2 = -e2 + B.eS + A.e4 
s3 = el + A,e3 - B,e4 
s4 = e2 + B,e3 + A.e4 

17. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 
calcul'^tion block comprises: 

2^-1 internal butterflies and a peripheral 
butterfly for the even values of index p, and 
2P-1 internal butterflies, otherwise. 

18. The calculation method according to claim 17, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the first 
transformation step with ]x being even, and otherwise 0. 

19. The calculation method according to claim 18, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0 and the 
seconda'ry mode signal is 1: 



s3 = e4 - e3 
s4 = e3 + e4 

2) if primary mode/ secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the' primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

si = (e3 + e4) - (el + e2) 

s2 = el " e2 

s3 = e4 - e3 

s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A,e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1 : 

si = el - A,e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 - B.e4 

s4 = e2 + B.e3 + A.e4 
20. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 



four inputs for receiving input samples and 
four outputs for providing output samples, 

- 'four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
5 inputs, 

in order to selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 

10 and a coefficient admitted on said corresponding 
additional inputs, 

and in that the final step furthermore performs an 
addition and subtraction between the first and the last 
output sample provided in the last transformation step. 

15 21. The calculation method according to claim 20, 

characterized in that, in transformation step Ep, each 
calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies. 

22. The calculation method according to claim 21, 

20 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
25 si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
30 permutation signal is 0: 

si = el + A.e3 - B.e4 



s2 = e2 + B.e3 + A.e4 
s3 == el - A. e3 + B. e4 
s4 = "e2 + B.eS + A.e4 
3) if the primary signal is 1 and the permutation 
5 signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 -B.e4 
s4 = e2 + B.e3 + A-e4 
10 23, The calculation method according to claim 9 or 

2 0, characterized in that the first and second binary 
addresses of p bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
15 address being greater than the first binary address. 

24. The calculation method according to claim 23, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

20 25. The calculation method according to claim 2 3 

or 24, characterized in that, if a peripheral butterfly 
is involved, the p+2 low-order bits of the first 
address are equal to 0, and the p+2 low-order bits of 
the second address form a number equal to 2^'^^-l. 

25 26. The calculation method according to claim 24 

or 25, characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly, be they peripheral or internal, are obtained 
by inverting the- (p+2) low-order bits of said first and 

30 second produced addresses. 
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27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 
odd-nuitibered address samples are stored in two separate 
memories. 

28. The calculation method according to claim 25, 
characterized in that the value of the parameter s of 
the coefficient assigned to an internal butterfly in 
transformation step Ep is coded by ia-2 bits, and is: 

if p+l=p-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, 

- if p+l<iJi-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, followed by y-p-S zero 
bits at the end of the number, 

if p+l>|i-2, the number formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, minus its p-p-1 low-order 
bits. 

29. The calculation method according to claim 4, 
in turn dependent on claim 3, in turn dependent on 
claim 2, characterized in that in each transformation 
step Ep, the butterflies are distributed among 2^ 
calculation blocks, 

in that each calculation block comprises one 
peripheral butterfly and internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank I'^'^a, 2^-Pa+2^'P-^-l, 2^"Pa-*-2^"P" 
^, 2^"^a+2^"^-l into output samples of the same rank. 
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and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
transforms the input samples of rank 2^"'^a+2x+l, 
Pa+2x+2, 2^-Pa-f-2^-P-^-.2x-3, 2^-Pa+2^-P-2x-2 into output 
samples of the sarae rank. 

30. The calculation method according to claim 29, 
characterized in further comprising a final step of 
modifying the ''sequence of the output samples provided 
in the last transformation step and classifying them in 
ascending order of index n. 

31. The calculation method according to claim 29 
or 30, characterized in that each butterfly is assigned 
a coefficient W^, whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
U2ns/u) ^^^^ g ^ [0_N/4-l] for a fast Fourier transform 
and is equal to e^'^ics/N) ^^^j^ ^ ^ [0.. N/4-1] for an 
inverse fast Fourier transform. 

32. Calculation method according to claim 31, 
characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 
is assigned coefficient with 5 = {x+l).2^. 

33. The calculation method according to claim 32, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

four inputs for receiving input samples and 
four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs. 



in order to . selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs. 

34. The calculation method according to claim 33, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the odd values. 

35. The calculation method according to claim 31 
or 34, characterized in that the secondary mode signal 
is 1 if the butterfly, be it peripheral or internal, is 
used for the first transformation step, and otherwise 
0. 

36. The calculation method according to claim 35, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode- and secondary mode signals 
are 0 : 

si = (el + e2) /2 
s2 = (el - e2)/2 
s3 = (e4 - e3) /2 
s4 = (e3 + e4) /2 

2) if the primacry mode signal is 0 and the 
secondary mode signal is 1: 



0 5C 0 



Si = [ {el + e4) /2-e2] /2 
s2 = [ (el + e4) /2-e2] /2 
s3 = [e3- (el-e4) /2] /2 



s4 = [e3+(el+e4)/2]/2 
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3) if the primary mode signal is 



1 



and 



the 



permutation signal is 0: 
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si = (el+e3) /2 

's2 = (e2+e4) /2 

s3 = (el-e3) .A - (e2 + e4) .B] J 2 

s4 = [-(el"e3).B + (e2+e4)^A]/2 



- 4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [(el-e3).A- (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 



37 • The calculation method according to claim 33, 
characterized in that the first and second binary 
addresses of yi bits are generated for each butterfly, 
2 0 each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 

38. The calculation method according to claim 37, 
characterized in that said first and second binary 

25 addresses are consecutive and an internal butterfly is 
involved. 

39. The calculation method according to claim 37 
or 38, characterized in that, if a peripheral butterfly 
-is involved, the |a-p low-order bits of the first 

30 address are equal to 0, and the p-p low-order bits of 
the second address form a number equal to N/2^^^-1. 
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s3 = (el+e3)/2 
s4 = (e2-e4) /2 



40. The calculation method according to claim 38 
or 39, characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly are obtained by inverting the \i-p low-order 
bits of both produced addresses. 

41. The calculation method according to claim 40, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories . 

42. The calculation method according to claim 41, 
characterized in that the value of the parameter s of 
the coefficient assigned to an internal butterfly in 
transformation step Ep is coded by \i-2 bits, and is: 

if |i-p-l=ia-2, the number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, 

if p-p-Kp-2, the number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, followed by p-1 zero bits at the 
end of the number, 

if u-p-l>vi-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, minus its p+1 low-order bits. 
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METHOD OF CALCULATING THE FAST FOURIER TRANSFORM AND 
THE INVERSE FAST FOURIER TRANSFORM 

This invention relates to a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of real numbers or _^a series of 
conjugated complex samples. 

The -Fourier transform is probably one of the most 
important tools for analyzing, designing and 
implementing signal processing algorithms, and the 
existence of efficient algorithms, such as that of the 
fast Fourier transform, has been a major factor for 
this situation. Although most Fourier transform 
algorithms are designed for transforming series of 
complex numbers, there are, however, various 
applications, such as image or acoustic signal 
15 processing or certain types of multicarrier modulation 
wherein the series to be transformed are real numbers. 

In general, the direct Fourier transform and the 
inverse Fourier transform respectively set up the 
following relations between two series of N complex 
20 numbers, x(n) and X (n) : 

N-l jjjl^^ 

X(n) = 2] J^(J<:)w''" with n e [0 ... N-l] and w"" = e'"*"^ 

lc=0 

2^ N-l 

= -Y.X{k)w-''" with n e [0 ... N-l] 

N )c=o 

In 1965, J.w. Cooley and J.W. Tukey described an 
algorithm allowing to rapidly calculate the Fourier 
25 transform of a series of N complex numbers in an 
article entitled "An Algorithm for the Machine 
Calculation of Complex Fourier Series", Math. 



computation. Vol. 19, 1965, pp. 297-301. This algorithm 
is very interesting when N is a power of two because it 
is then particularly straightforward to implement. This 
algorithm requires p calculation steps where u=log2 (N) . 
It is based on breaking down the series to be 
transformed into two interleaved sub-series. There are 
two kinds of interleaving: time interleaving and 
frequency interleaving. Both kinds of interleaving are 
explained more in detail in the course of the 
description with reference to Figs. lA and IB. 

Specific algorithms have been developed from this 
basic algorithm to deal with the case of real number 
series. The Fourier transformation of a series of 8 
real numbers according to a time interleaving algorithm 
and according to a frequency interleaving algorithm is 
illustrated in Figs. lA and IB. For each time 
interleaving Fourier transform algorithm there is a 
frequency interleaving algorithm, which corresponds to 
a double inversion of the series of transformation 
operations, on the one hand, and for each butterfly 
circuit, of the proper transformation operations, on 
the other hand. Whatever the interleaving chosen, the 
transformation method requires three transformation 
steps EO, El, and E2, these steps being implemented 
through a set of four butterfly circuits CC, generally 
known as "butterflies" in technical speak. Each 
butterfly, represented in the figure by a point of 
intersection between two columns of numbers, performs 
calculations on two numbers, real or complex. The 
symbols R and C respectively identify a real number and 
a complex number. The sequence of real and/or complex 



numbers after the transformation steps depends on the 
interleaving chosen. 

The time interleaving algorithm is generally 
chosen for calculating the Fourier transform of a 
series of real numbers because of the symmetrical 
distribution of . real and ..complex- numbers throughout the 
steps. On the other hand, the frequency interleaving 
algorithm is more suitable for the direct or inverse 
transformation of a series of conjugated complex 
numbers . 

When the series to be transformed x(n) is real/ 
the Fourier transform verifies the following relation: 

x(n) is real if and only if 

X(n) = X* (-n)=X* (N-n) ; (D 
where * designates the conjugating operation. 

For a series x(n) of N real nximbers, the following 
results are inferred from this relation: 

- X(0) and X (N/2) are real; 

- X(n) = X* (N-n) for 1 ^ n ^ N/2 -1 

Relation (1) highlights the presence of redundant 
information in the X(n) series. 

It should be noted that the transformation method 
is generally implemented by a single set of 
butterflies, the operating mode of which is modified as 
the transformation goes along. At each change of 
operating mode, results are stored in a memory having N 
storage locations, the output samples of a butterfly 
replacing the corresponding input samples of the same 
rank in .the memory. This method of applying the 
algorithm is generally known as the "in place" method. 
This method has a major advantage: if the elements of 
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the x(n) series are processed in the first 
transformation step in bit-reversed order of index n, 
the numbers of the X(n) series are output in the last 
transformation step in ascending order of index n and 
vice versa, 

A known transformation method is shown by way of 
example in Fig, 2. This method performs the Fourier 
transformation of a real x (n) series according to a 
complex time interleaving algorithm. In this example, 
the x{n) series to be transformed comprises sixteen 
real samples, x{0) to x(15). The transformation method 
comprises four transformation steps Ep with 0^p<3. The 
samples of the x(n) series are shown in the first 
transformation step in bit-reversed order of their 
15 index n. 

At this stage of the explanations, the terms used 
in the course of the description should be defined. The 
rank of a sample is taken from the position it occupies 
in the series of samples to which it belongs* The index 
20 of a sample then corresponds to the starting rank of 
this sample. 

The intermediate results obtained in the various 
transformation steps are represented by the series 
A(n), B(n), and C(n). The samples of the series x(n), 
25 A(n), B(n), C (n) , and X(n) are stored in double storage 
locations, one storage location being reserved for the 
real portion of the sample and the other location being 
reserved for the imaginary portion thereof. AR(n) and 
Ai (n) respectively designate the real portion and the 
imaginary portion of the index n sample of the A(n) 
series. Butterflies are represented in the figure by 



points of intersection between columns of storage 
locations. Each butterfly is assigned a coefficient 
symbolized in Fig. 2 by a pair of coordinates A/B where 
A and B respectively designate the real portion and the 
imaginary portion of the coefficient W^. Coordinates 
1/0 and 0/-I are respectively assigned to coefficients 
W°=l and w^/''=W^=-j. For the sake of clarity and in 
order to simplify their formulation, the remaining 
coefficients have been represented by the following 
pairs : 

-> 2 / -4 -> -4 / -2 

-> 3 / -3 -> -3 / -3 

-> 4 / -2 W'' -> -2 / -4 

These pairs of coordinates are graphically 
represented in Fig. 3. In fact, coordinates A and B 
respectively represent a cosine value and a sine value. 
This coefficient participates in the calculation 

performed by the butterfly. Furthermore, the 
butterflies are distributed at each transformation step 
among N/2P*^ calculation blocks, each calculation block 
comprising 2^ butterflies. In the course of the 
description, the parameter q designates the rank of the 
calculation blocks within the same transformation step; 
q is included between 0 and (N/2P*^)-1. 

In the first transformation step Eo/ the 
butterflies are distributed among eight calculation 
blocks, each comprising a butterfly performing an 
operation. on two complex or real samples. If el and e2 
are to designate the samples applied to the inputs of a 
butterfly, the latter outputs samples si and S2 defined 
as follows: 
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Si = el + W^.e2 and s2 = el -(W^.e2) 
where is the coefficient assigned to said butterfly. 
For this first transformation step, coefficient 
W*^=l is assigned to the eight butterflies. As the 
5 samples x(n) and the coefficient W° are real, the 
samples A(n) obtained at the end of step Eo are real* 

For the second transformation step, Ei, the 
butterflies are distributed among four calculation 
blocks each "comprising two butterflies. Coefficient 

10 W*^=l is assigned to the first one of these butterflies; 
thuS/ the first butterfly of each calculation block 
provides two real samples. The second butterfly of the 
calculation blocks is associated with coefficient 
W^^^=W*=-j and generates two conjugated complex samples. 

15 The output samples obtained at the end of step Ei. are 
designated by the series B(n}. 

For the third transformation step, Ez, the 
butterflies are distributed among two calculation 
blocks each comprising four butterflies respectively 

20 associated with coefficients W^, and W^. The 

output samples in step E2 are designated by the series 
C(n), Finally, for the fourth transformation step, E3, 
a single calculation block comprising eight butterflies 
respectively associated with coefficients W°, W^, W^, 

25 W^, WS W^, W^, and w"', is provided. This transformation 
step generates the transformed series X (n) . 

Given relation (1), the X(n) series comprises, on 
the one hand, real samples, X(0) and X(8), and on the 
other hand, complex samples, X(l) to X(7) and X(9) to 

30 X(15), samples X(15) to X(9) respectively being the 
conjugates of samples X(l) to X(7). The X(n) series 



therefore contains redundant information. The storage 
locations outlined in thick stroke in Fig. 2 designate 
the storage locations enclosing the conjugated values 
of the complex samples contained in the storage 
5 locations associated therewith by an arrow. The 
intermediate results series B(n) and C(n) also contain 
redundant information . 

It is then possible to delete this redundant 
information in order to reduce by half the size of the 
10 sample storage memory as well as the number of 
butterflies. 

However, removing redundant information stored in 
the storage locations outline in bold strokes in Fig. 2 
implies the complete reorganization of the 

15 transformation steps of Fig. 2. Thus reorganizing the 
transformation has the effect of modifying the output 
sequence of samples X{n). 

The problem therefore consists in reducing the 
size of the storage memory and the number of 

20 butterflies while maintaining the output sequence of 
the X(n) samples. It is the object of the invention to 
offer a method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
series of N real samples x(n), with N power of 2, 

25 operating according to a time interleaving algorithm, 
which provides the series of samples X(n) in ascending 
order of index n and uses limited calculation and 
storage means. 

For this purpose, the object of the invention is a 

30 method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 




defined by a series of N real starting samples x(n), 
with N power of two and n e [0,.N-1], comprising 
successive transformation steps for transforming input 
samples into output samples, all the transformation 
5 steps being performed by means of a single set of 
butterflies with several inputs and several outputs^ 
the operating mode of which is modified selectively in 
each transformation step, the input and output samples 
of each transformation step being stored in a storage 
10 memory, a series of N output samples y(n) 
representative of the fast Fourier transform or the 
inverse fast Fourier transform of the starting samples 
x(n) being provided in the last transformation step, 

characterized in that output samples y(n) are 

15 real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x (n) processed in the first transformation step 
20 are classified in bit-reversed order of their index n, 
output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 
25 y(0) - Re[X(0) ] 

y(n) = Re [X { (n+1) /2) ] for n being odd and 

different from N-1 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

30 y (N-1) = Re [X (N/2) ] 




where the X(n) samples, with n g [0..N-1] designate 
the complex samples of the series corresponding to the 
fast or inverse fast Fourier transform of the series of 
starting samples x(n), 
5 For methods operating according to a frequency 

interleaving algorithm, the invention also relates to a 
method of calculating the fast Fourier transform or the 
inverse fast Fourier transform of a digital signal 
defined by a series of N complex samples X(n) 
10 conjugated by pairs represented by a series of N real 
starting samples y(n), with N power of two and n e 
[0..N-1], the starting samples y(n) being defined as 
follows: 

y(0) = Re[X(0) ] 

15 y(n) = Re[X( (n+l)/2) ] for n being odd and 

different from N-1 
y(n) = Im[X{n/2)] for n being even , and 

different from 0 

y(N-l) = Re[X(N/2) ] 

20 this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of N output samples x(n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 

25 step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
is modified selectively in each transformation step, 
the input and output samples of each transformation 

30 step being stored in a storage memory. 
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characterized in that output samples x(n) are 

real/ 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output -samples x(n) are output in the last 
transformation step in bit-reversed order of index n. 

The inventive calculation methods perform 
operations on real samples and consequently ..use limited 
calculation and storage means in comparison with the 
method shown in Fig. 2. 

According to another feature of the invention, in 
each transformation step, the butterflies transform 
input sample pairs, the ranks of the input samples of 
the same pair within the series of input samples of 
said transformation step being symmetrical with respect 
to a center between the end rank values of the input 
samples transformed by said butterfly. Input samples 
processed by the same butterfly are thus symmetrically 
linked together by pairs. The result is simplified 
handling of sample addressing. 

According to another aspect of the invention, the 
method preferably comprises p-1 transformation steps Ep 
with ia=log2(N) and p e [0..U-2]. 

Other features and advantages of the invention 
will be apparent from reading the following detailed 
description, which is made with reference to the 
appended drawings, where: 



Figs. lA and IB, already described/ 
respectively represent a Fourier transformation of 
eight real numbers according to a time interleaving 
algorithm and according to a frequency interleaving 
5 algorithm; 

Fig. 2/ already described, illustrates the 
transformation of a series of 16 real numbers into a 
series of 16 complex numbers according to a complex 
time interleaving algorithm; 
10 - Fig. 3, already described, graphically 

represents the . mapping of coefficients and 
coordinate pairs A/B; 

Fig. 4 represents modifications applied to part 
of the transformation of Fig. 2; 
15 - Fig. 5 illustrates a modified transformation 

only processing real nximbers; 

Fig. 6 is a representation of the method of 
calculating the fast Fourier transform according to the 
invention; 

20 - Figs, 7A and 7B respectively illustrate 

permutations performed on peripheral butterflies and on 
internal butterflies with odd-numbered rank of the 
transformation of Fig. 5; 

Fig. 8 represents an embodiment of the 

25 transformation method according to the invention, 
comprising ]i transformation steps; 

Fig. 9 illustrates a grouping of peripheral 
butterflies according to a first embodiment of a 
transformation, method comprising p-1 transformation 

30 steps; 
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Fig. 10 represents a first embodiment of a 
transformation method comprising p-1 transformation 
steps; 

Fig. 11 represents a butterfly design relating 
to the transformation method illustrated in Fig. 10; 

Fig. 12 represents an alternative of the 
embodiment of Fig. 10; 

Fig. 13 represents a butterfly design relating 
to the embodiment of Fig. 12; 

Fig. 14 illustrates a grouping of peripheral 
butterflies according to a second embodiment of a 
transformation method comprising ia-1 transformation 
steps; 

- Fig. 15 represents a second embodiment of a 
transformation method comprising p-1 transformation 
steps, with p being even; 

Fig. 16 represents an alternative of the 
preceding embodiment, with p being odd; 

Fig. 17 represents a butterfly design relating 
to the embodiments shown in Figs. 15 and 16; 

Fig. 18 represents a third embodiment of a 
transformation method comprising p-1 transformation 
steps; 

Fig. 19 represents a butterfly design relating 
to the embodiment of Fig. 18; 

Fig. 20 represents the addresses that are 
associated with the various butterflies implemented in 
the embodiment shown in Fig. 12; 

~ Fig. 21 represents the addresses that are 
associated with part of the samples of a transformation 
method processing a series of 32 real samples; 
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Fig. 22 represents a sample embodiment of a 
transformation method operating according to a 
frequency interleaving algorithm; 

Fig. 23 represents a butterfly design relating 
to the embodiment of Fig. 22; 

Fig. 24 represents the addresses that are 
associated with the various butterflies of the 
embodiment illustrated in Fig. 22. 

According to the invention, only one part of the 
X(n) samples is calculated, the other part of the 
samples being redundant. E.g., the calculation could be 
limited to the X (n) numbers for 0<n<N/2. However, a 
more adequate solution is calculating the X (n) niombers 
for 0<n<N/4 and N/2<n<3N/4. The latter solution is 
preferable because it only involves the first (N/4)+l 
butterflies for calculating X(n), as the last (N/4)-l 
butterflies can be removed- The storage locations thus 
released can be used for storing the real portion or 
the imaginary portion of the remaining X (n) numbers . 
The size of the storage memory can thus be halved if 
the size of the storage locations is limited to storing 
a real number instead of a complex number. The real 
portion of the complex X(n) number is stored in the 
storage location that has been assigned thereto 
initially whereas the imaginary portion is stored in 
the storage ' location initially assigned to the number 
X(N-n) . The same operation can be performed for the 
series of intermediate results B(n) and C(n). The 
transformation method of the invention is limited to 
calculating N real samples y(n) instead of N complex 
samples. The real samples y{n) are defined as follows: 



y(0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-1 

y(n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = Re[X(N/2) ] . 

Removing redundant information and reorganizing 
storage - locations deeply modifies the diagram of 
Fig- 2* Fig/ 4 illustrates the reorganizations 
performed on the part relating to the calculation of 
samples A(l) , A(3) , A(5) , A(7), B(l), B(3), B(5), B(7), 
C(l)/ C(3), C(5), and C(7). Redundant samples to be 
removed are B(3), B(7), C(5), and C(7). When the size 
of the storage locations has been reduced and the 
samples remaining in theses storage locations have been 
reorganized, the butterflies perform calculations on 
real numbers. The butterflies to which the pair 0/-1 is 
assigned perform calculations on two real numbers. In 
practice, they copy onto their first output the number 
present at their first input and multiply by -1 the 
number present at their second input and provide it to 
their second output. The butterflies to which the pair 
1/0 is assigned perform an addition and a subtraction 
on two real numbers. Finally, the other butterflies 
perform operations on four real numbers. 

The transformation thus rearranged is illustrated 
in Fig. 5. In this figure, the butterflies associated 
with the pairs 1/0 and 0/-1 corresponding to the 
coef f idients and are connected with storage 

locations by thick strokes. This figure shows that the 
reorganization of the transformation method steps 




modifies the output sequence of the y(n) samples and 
therefore the output sequence of the X(n) series. 
Furthermore, this transformation method no longer has a 
specific symmetry allowing to link by pairs the samples 
5 processed by the same butterfly. The result is a very 
complicated address management of the samples to be 
applied to the butterfly inputs. 

Fig.' 6 shows the steps of a method of calculating 
the fast Fourier transform or the inverse fast Fourier 

10 transform of a series of N real numbers x(n), with N 
power of 2, operating according to a time interleaving 
algorithm. It mainly comprises transformation steps 2 
for transforming N starting samples x(n) classified in 
bit-reversed order of their index n into real output 

15 samples y(n) representative of this Fourier transform 
classified in ascending order of index n. 

Advantageously, it comprises a preliminary step 1 
for ranking the N real starting samples x(n) to be 
transformed in the bit-reversed order of their index n 

20 if the samples x(n) are not already in this sequence 
and a final step 3 for generating the N complex samples 
X(n) corresponding to the fast Fourier transform of the 
starting samples x(n) from the N real samples y(n) 
obtained at the end of the transformation steps. 

25 The methods that will be detailed in the course of 

the description will be more in particular for 
calculating the fast Fourier transform of a real 
series. Also, the coefficients assigned to the 

butterflies for implementing the inventive method will 

30 be of the type e"^*^''^''^' with s>0. For calculating the 
inverse fast Fourier transform the calculation method 




is the same, however, the coefficient is of the type 

In order to obtain at the same time y(n) samples 
sorted in ascending order of index n and symmetry of 
5 calculation, according to the invention, it is 
suggested to modify the calculations performed by the 
butterflies of the odd rank q calculation blocks of the 
transformation illustrated in Fig. 5 in accordance with 
the diagrams of Figs. 7A and 7B, 

10 As the butterflies associated with coefficient 1/0 

(Fig* 7A) of the odd rank calculation blocks are 
involved, provision is made for permutation both 
outputs of the butterfly and multiplying by -1 the 
result provided at the second butterfly output. 

15 As the four input butterflies (Fig. 7B) are 

involved, provision is made for permutation the first 
two outputs with the last two ones. 

This method is applied to the whole transformation 
and then a method is obtained providing at the output 

20 y(n) samples in ascending order of index n. This method 
is illustrated in Fig. 8. The symbol^ placed above the 
calculation blocks designates the calculation blocks 
wherein butterflies have been modified, i.e. odd rank 
calculation blocks. Due to the symmetry of calculation 

25 of the transformation method, the intersecting points 
representing the butterflies are superposed inside each 
calculation block. 

The butterflies to which the coefficient 1/0 is 
assigned are called peripheral butterflies because they 

30 perform calculations on samples arranged at the ends of 
the calculation block. The other butterflies are called 




internal butterflies. It should be noted that in each 
transformation step, not all the samples are processed 
always and that unprocessed samples are kept in their 
storage locations to be processed in subsequent steps, 
5 or else produced as output if they are already in their 

final shape . . 

The transformation method thus modified provides 
y(n) samples in ascending order of index n and, in each 
transformation step, has symmetry of calculation 
10 facilitating the . addressing of the samples to be 
processed. 

According to another aspect of the invention, the 
calculation method advantageously comprises p-l 
transformation steps. Several embodiments derived from 
15 that of Fig. 8 and comprising ii-l transformation steps 
are therefore shown in the course of the description. 
One butterfly design is associated with each of these 
embodiments. 

All of these embodiments have the following 
20 features in common: 

in each transformation step, provision is made 
for calculation blocks and each calculation block 

comprises an peripheral butterfly and/or 2^-1 internal 
butterflies; all butterflies, be they peripheral or 
25 internal, perform calculations on four real samples; 

- the ranks of the samples processed by the same 
butterfly are defined as follows: if, in transformation 
step Ep, a peripheral butterfly belonging to the rank a 
calculation block is taken into consideration, it 
30 transforms the input samples of rank 2P^^a, 2^*^a+2P'"^-l, 
2^^^a-^2^''^ , 2^''^a+2^''^-l into output samples of the same 




rank, and, if an internal rank x butterfly in a rank a 
calculation block is taken into consideration in step 
Ep, it transforms the input samples of rank 2P*W2t+1, 
2P^^a+2x^-2, 2P*W2P^^-2x-3, 2P^^a-.2P^^-2x-2 into output 
samples of the same rank, with P>1; 

- the coefficient assigned to the internal rank x 
butterfly of the rank a calculation block in step Ep is 
equal to with 6 = (x+1) , . 

In all of these embodiments, the input samples of 
each butterfly form sample pairs, the sample ranks of 
the same pair in the series of input samples of a 
transformation step being symmetrical with respect to 
the center value of the end ranks of the input samples 
transformed by said butterfly. This center corresponds 
to the value 2P^2^+2P^^-l/2 . Therefore, one just has to 
know the rank of two of the four samples to be applied 
to the inputs of the butterfly in order to infer 
therefrom the rank of the other two. Addressing these 
samples is thus simplified. This will be explained more 
in detail at a later point of the description. 

Thus, according to a first embodiment, the 
neighboring calculation blocks in each transformation 
step are grouped by pairs. The peripheral butterflies 
of the same calculation block are then merged into a 
single peripheral butterfly. A sample merging of two 
peripheral butterflies is shown in Fig. 9. This example 
relates to peripheral butterflies associated with 
samples x(0), x(8), x(4), andx(12). 

Furthermore, as the transformation steps do not 
process all samples each time, certain calculations can 



be anticipated. E.g., calculating samples Cr{1)/ Ci(1), 
Cr(3), and Ci(3) can be done in the second 
transformation step. The result is the diagram of Fig. 
10 showing a first embodiment of the transformation 
circuit wherein the method only comprises y^-l 
transformation steps. In transformation step Ep, each 
calculation block, has a peripheral butterfly and 2^-1 
internal - butterflies. It can be considered that this 
method only ^comprises 3 transformation steps, the 

fourth step being limited to performing an addition. and 
a subtraction. This addition and this subtraction are 
preferably performed during the final step, and in 
order to limit the number of transformation steps. 

A butterfly design associated with the embodiment 

of Fig. 10 is represented in Fig. 11. It comprises: 

four inputs for receiving input samples el/ e2/ 

eS/ e4, and four outputs for providing output samples 

si, s2, s3, s4, and 

- three additional, respectively primary mode MP, 

permutation PERM, and coefficient COEF, inputs. 

This butterfly is responsible for selectively 

applying to input samples el, e2, e3, and e4, various 

transformation operations each determined by the values 

assigned to primary mode, permutation signals and to 

coefficient admitted at the corresponding additional 

inputs . 

The primary mode signal is 0 for a peripheral 
butterfly and 1 for an internal butterfly. When the 
permutation signal is a 1, the output samples si and s2 
of the butterfly are swapped with output samples s3 and 
s4. This permutation is only possible if the butterfly 



is an internal one. Finally, the coefficient 
associated with the butterfly is applied to the 
coefficient input COEF, 

Thus, if the complex coefficient W^=A+j.B is 
5 applied to the coefficient entry of the butterfly, the 
latter provides the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
si = el + e2 

10 s2 = el - e2 

s3 = e4 - e3 
. s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
permutation signal is 0: 

15 si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B,e3 + A.e4 

3) if the primary signal is 1 and the permutation 
20 signal is 1: 

si = el - A.e3 + B,e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el + A.e3 -B,e4 

s4 = e2 + B.e3 + A,e4 
25 According to an alternative embodiment, provision 

can be made for addition and subtraction to be 
performed inside the peripheral butterfly of the last 
transformation step. This embodiment is shown in Fig. 
12. For ' this purpose, the corresponding butterfly 
30 design has a fourth additional input called secondary 
mode input MS to which a secondary mode signal is 



applied. This signal is 1 for the peripheral butterfly 
of the last transformation step, otherwise it is 0. 
This design is illustrated in Fig. 13. This design has 
an additional operating mode in comparison with the 
5 preceding one; thus, when the primary mode signal is 0 
and the secondary mode signal is 1, the output obtained 
is : 

si = el + e2 + e3 + e4 
s2 = el - e2 
10 s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

According to a second embodiment derived from the 
diagram of Fig. 8, the neighboring peripheral 
butterflies in even index p transformation steps Ep are 

15 grouped by pairs and are merged with the peripheral 
butterfly of the second odd index step in order to form 
a new peripheral butterfly at the odd index step. This 
grouping is illustrated in Fig. 14 in an example. In 
this example, the peripheral butterflies of the first 

20 transformation step processing samples x(0), x(8), 
x(4), and x(12) are merged with the peripheral 
butterfly of the second step processing samples Ar(0) 
and Ar(2) . Also, the peripheral butterflies of the 
first transformation step processing samples x{2), 

25 x(10), x(6), and x(14) are merged with the peripheral 
butterfly of the second step processing samples Ar(4) 
and Ar(6). The two butterflies obtained are different 
in that the second one performs in addition a 
permutcctibn between the first and second outputs. If 

30 this grouping is applied to the whole transformation 
illustrated in Fig, 8, the result is that step Eo no 



longer uses any butterflies and can be removed. The 
resulting transformation method is illustrated in Fig. 
15. 

However, two cases should be distinguished for 
5 this transformation method: the case where N is an even 
power of two (u being even) and the case where N is an 
odd power of two {]x being odd) . 

In case y is even, there is an even number of 
transformation steps in the embodiment shown in Fig. 8 
10 and grouping the peripheral butterflies of even index 
steps with those of the following odd index steps is no 
problem. This case corresponds to the diagram in Fig. 
15. 

In case ]x is odd, the peripheral butterfly of the 

15 last even index step cannot be grouped with other 
peripheral butterflies. Therefore, a specif ic . operating 
mode should be provided for this case. This case is 
illustrated in Fig. 16, this figure representing the 
transformation of a series of eight real samples (ii=3) . 

20 The peripheral butterfly of the last transformation 
step of this circuit could not be merged with other 
peripheral butterflies. 

The butterfly design associated with this second 
embodiment is illustrated in Fig. 17/ it differs from 

25 the preceding design in that the secondary mode signal 
is 1 when a peripheral butterfly for implementing the 
last step is involved and ]i is odd, and in that 
permutation applies to all the butterflies of the even 
rank calculation blocks, 

30 The calculations performed by the butterfly are 

also different and are defined as follows: 



1) if primary mode/ secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
5 s3 = e4 - e3 

s4 = (el + e2) - (eS + e4) 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e4 
10 s2 = e2 

s3 = e3 
s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

15 si = (e3 + e4) - (el + e2) 

s2 = el - e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
20 permutation signal is 0: 

si = el + A,e3 - B.e4 

s2 = e2 + B.e3 + A.e4 

s3 = el - A.e3 + B.e4 

s4 = -e2 + B,e3 + A,e4 
25 5) if the primary mode signal is 1 and the 

permutation signal is 1: 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
/ s3 = el + A.e3 - B.e4 
30 s4 = e2 + B.e3 + A.e4 
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In case p is odd, it is also possible on the one 
hand to provide for the peripheral butterflies to be 
grouped for implementing the first transformation step 
in the same way as in the first embodiment, and on the 
5 other hand, for the butterflies of the other steps to 
be grouped as in the third embodiment. Grouping 
butterflies from the second step on is then performed 
by assigning an even index to the first transformation 
step. These groupings are represented in Fig, 18. 

10 The butterfly design corresponding to this 

embodiment is represented in Fig. 19. The secondary 
mode signal is 1 for a peripheral butterfly 
implementing the first transformation step of the 
circuit and if ]i is even. The calculations performed by 

15 this butterfly are the following ones: 

1) if the primary mode signal is 0 and the 
secondary mode signal is 1 : 

si = el + e2 
s2 = el - e2 
20 s3 = e4 - e3 

s4 = e3 + e4 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
25 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

30 si = (e3 + e4) - (el + e2) 




s3 = e4 - e3 

s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

5 si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A,e4 
s3 = el - A.e3 -f B.e4 
s4 = -e2 + B,e3 + A.e4 

5) if th'e primary mode signal is 1 and the 
10 permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 
15 In all the embodiments described before, the ranks 

of samples processed by the same butterfly are 
symmetrical by pairs with respect to a center value. 
One just has to know the rank of the first two input 
samples of the butterfly to infer the other two ones 
20 therefrom by syrmnetry. If the input and output samples 
are saved in storage locations the address of which 
corresponds to the rank of these samples, addressing 
the latter will be simplified. 

Indeed, all that needs to be done then is to 
25 generate two addresses per butterfly, as the other two 
ones can be inferred by symmetry. Furthermore, it 
should be noted that the addresses of the input samples 
and those of the corresponding output samples are the 
same because the transformation is applied according to 
30 an in place method. 



The addresses associated with the various 
butterflies relating to the transformation method of 
Fig. 12 are shown in Fig. 20. The address of a sample 
is taken to be equal to the rank of this sample in the 
sample series to which it belongs. The series to be 
transformed in the example of Fig. 12 comprises 16 
samples. Therefore, 16 addresses have to be produced, 
with four-bit addressing. For a series of N samples, 
log2{N)-bit addressing is performed. 

Each transformation step of the method of Fig. 12 
is implemented. by four butterflies each processing four 
real input samples. The binary addresses of the samples 
to be processed in each step are therefore distributed 
among four groups of four addresses. The address groups 
relating to one peripheral butterfly are contained in 
bold line boxes and the address groups relating to an 
internal butterfly are contained in thin line boxes. 
Furthermore, the address groups relating to the same 
calculation block are grouped in dotted line boxes. 

The first peripheral butterfly in step Eo 
processes the first four samples of the x(n) series 
resulting from the preliminary classification step. The 
corresponding addresses to be generated for this 
butterfly are therefore 0000, 0001, 0010, and 0011. 
Also, considering the first peripheral butterfly of 
step El, it processes the samples of rank 0, 3, 4, and 
7 of the A(n) series. The corresponding addresses to be 
generated for this butterfly are therefore 0000, 0011, 
0100, and 0111. 

The (ia-2-p) high-order bits of the addresses 
indicate rank q of the calculation block to which the 



address is related. Thus, the two high-order bits of 
the addresses generated for the first calculation block 
of step Eo are 00. It should also be noted that, as the 
last step only uses a single calculation block/ the 
5 addresses generated for this calculation . block have no 
bit relating to the rank of this block (u-p-2-0) . 

For each butterfly, provision is made for 
generating only two binary addresses; the other two 
ones are obtained by inverting the- (p+2) low-order bits 

10 of the generated addresses* E.g., considering the 
addresses of the first peripheral butterfly in step Ei, 
only addresses 0000 and 0011 are generated and 
addresses 0111 and 0100 are obtained by inverting the 3 
low-order bits of the generated addresses. 

15 A first and a second address per butterfly are 

then produced, these addresses being consecutive for an 
internal butterfly. When a peripheral butterfly is 
involved, the p+2 low-order bits of the first address 
are equal to 0,. and the p+2 low-order bits of the 

20 second address form a number equal to 2^''^-l. 

Regarding the coefficients to be applied to the 
coefficient inputs COEF of the internal butterflies, 
they are stored in a memory of the calculation circuit. 
Only N/4 coefficient values are required for 

25 calculating the Fourier transform, ia-2-bit addressing 
is performed for these coefficients. In the embodiments 
described previously where N=16, only the pairs 1/0, 
2/-4, 3/-3, and 4/-2 are used which correspond to 
coefficients W^, W^, W^, 

30 According to the invention, the address of these 

four coefficients must therefore be known. The address 



associated with each coefficient is chosen to be 

equal to the value of the power s. Consequently, the 
addresses of coefficients W^, W^, are 00, 01, 10, 

11/ respectively. 

So as not to be obliged to generate these 
addresses, according to the invention, the addresses 
produced for addressing the samples are used. The 
address of the coefficient which is assigned to a 
butterfly is "included in the second address produced 
corresponding to the greater one of both addresses • 

However, we must distinguish between three cases: 

a) when p+l=ii-2, the coefficient address 
corresponds to the number formed by the p+1 low-order 
bits of the second address generated for this internal 
butterfly. This is the case for the second step (p=l) 
in the example of Fig. 20. The 2 low-order bits of the 
second address are 10 and therefore designate 
coefficient W^, 

b) if p+l>)a-2, the coefficient address corresponds 
to the number formed by the p+1 low-order bits of the 
second address generated for this internal butterfly, 
minus its u-p-1 low^order bits. This is the case for 
the third step (p=2) in the example of Fig. 20. The 3 
low-order bits of the second address generated for the 
first internal butterfly are 010. When the last bit (ii- 
p-l=l) of this number is taken away, the number 01 is 
obtained relating to coefficient W^, This case always 
corresponds to the last step of a transformation method 
compris'ing ]x-l transformation steps. 

c) if p+Kyi-2, the coefficient address corresponds 
to the number formed by the p+1 low-order bits of the 
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second address generated for this internal butterfly, 
followed by u-p-3 zero bits at the end of the number. 
This case is illustrated in Fig, 21. This figure 
represents the addresses relating to an internal 
5 butterfly for implementing the second transformation 
step (p=l) of a transformation method designed for 
processing a series of 32 real samples (N=32 and p=5) . 
This. butterfly transforms the samples of rank 1, 2, 5, 
and 6 of the series of samples obtained at the end of 

10 the first step. The two low-order bits of the second 
generated address are 10 and if a zero (p-p-3=l) is 
added at the end of the niomber, the number 100 is 
obtained designating coefficient , 

ThuS/ both addresses produced by the address 

15 generator for one butterfly is used to address both 
four samples to be processed and the coefficient 
relating to the butterfly. 

Preferably, even address and odd address samples 
will be stored in two separate memories. Thus, it will 

20 be possible to read two input samples simultaneously, 
and it will be possible to write the resulting output 
samples simultaneously, which means saving process time 
for the series to be transformed. 

As mentioned before, for each calculation method 

25 operating according to a time interleaving algorithm 
there is a corresponding method operating according to 
a frequency interleaving algorithm. All that has to be 
done to obtain it is, on the one hand, to invert the 
series of transformation operations of the 

30 corresponding time interleaving method, and on the 
other hand, for each butterfly, to invert the 




transformation operations themselves with respect to 
those of the corresponding method. 

In addition, the invention also relates to a 
method of calculating the fast Fourier transform or the 
5 inverse fast Fourier transform of a series of N complex 
samples X (n) conjugated , by pairs, with- N power of 2, 
operating according to a frequency interleaving 
algorithm. The series of N complex samples X (n) is 
represented by a series of N real samples y(n) defined 
10 as follows: 

y (0) = Re[X(0) ] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-1 
y(n) = Im[X(n/2)] for n being even and 

15 different from 0 

y(N-l) = Re [X (N/2) ] 

According to the invention, this method 
substantially comprises transformation steps for 
transforming input samples into output samples- Real 

20 samples y(n) are processed in a first transformation 
step and the last step provides a series of N output 
samples x(n) representative of the fast or the inverse 
fast Fourier transform of the sample series X (n) . As 
for the time interleaving circuits, each transformation 

25 step is implemented by a set of butterflies with 
several inputs and several outputs. As the 
transformation is done according to an in place method, 
all the steps are performed by means of a single set of 
butterflies, the operating mode of which is modified in 

30 each transformation step. In each transformation step, 
input and output samples are stored in a storage 
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memory. When they have been transformed, the output 
samples of the same butterfly replace the corresponding 
input samples of the same rank in the storage memory. 

According to the invention, if the samples y(n) 
5 introduced in the first transformation step are 
classified in ascending order of index n, the output 
samples x{n) are ' provided in the last transformation 
step in bit-reversed order of index n. The output 
sequence of the x(n) samples can then be modified by a 
10 final step so as to classify them in ascending order of 
index n. 

An embodiment of such a transformation is 
represented in Fig. 22. It is inferred from the 
embodiment of Fig. 12 by inverting the functional 

15 arrangement of the transformation steps of Fig. 12 
(mirror image of what it is for time interleaving) . 
This embodiment allows the calculation of the inverse 
fast Fourier transform of a series of 16 real samples 
y(n) representative of a series of 16 complex samples 

20 X(n) conjugated by pairs. Coefficients are therefore 



This embodiment comprises three transformation 
steps Ep with 0<p<2 . In each transformation step Ep, 
the butterflies are henceforth distributed among 2^ 

25 calculation blocks, these calculation blocks being 
sorted in each step according to an ascending rank q 
from 0 to 2^-1. Each calculation block has a peripheral 
butterfly and N/2P^^-1 internal butterflies. The y(n) 
samples are applied in the sequence of index n in the 

30 first step. 



of the type e 
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The ranks of the samples processed by the same 
butterfly are defined as follows: in step Ep, 
considering a peripheral butterfly belonging to the 
rank a calculation block, it transforms the input 
5 samples of rank a^'^a, 2^-Pa+2^-P-'-l , 2^-Pa+2^-P-S 2^-Pa+2^- 
^-1 into output samples of the same rank, and, in step 
Ep, considering an internal rank x butterfly in the 
rank a calculation block, it transforms the input 
samples of rank 2^"Pa+2T+l, 2^"Pa+2x+2, 2^-Pa+2^"P-2x-3, 2^' 

10 Pa+2^"'^-2T-2 into output samples of the same rank. 
Finally, the coefficient assigned to the rank t 
internal butterfly of the rank a calculation block in 
step Ep is equal to with 5 = (x+1) .2^, 

The real samples x(n) obtained at the end of the 

15 method are provided in bit-reversed order of index n. 

At the butterflies, the coefficients are of the 
type e^*^"^^^* and inputs and outputs have been inverted 
with respect to the embodiment of Fig, 12. 
Consequently, the operations performed by the 

20 butterflies of this embodiment are different from those 
performed by the butterflies of Fig, 12. A butterfly 
design associated with the embodiment of Fig, 22 is 
represented in Fig. 23, 

Just like the design associated with Fig. 12, it 

25 has four data inputs and four data outputs as well as 
four additional inputs, respectively primary mode MP, 
secondary mode MS, permutation PERM, and coefficient 
COEF inputs. The primary mode signal is 0 for an 
peripheral butterfly and 1 for an internal butterfly, 

30 The permutation signal is 0 for even values of rank q 
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and 1 for odd values* Finally/ the secondary mode 
signal is 1 if the peripheral butterfly is used for 
implementing the first step, and otherwise 0. 

The calculations performed by this butterfly are 
the following ones (W^=A+j.B): 

1) if the primary mode and secondary mode signals 
are 0: 

si = (el + e2) /2 
s2 = (el - e2) /2 
s3 = (e4 - e3) /2 
s4 = (e3 + e4)/2 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = [ (el+e4) /2+e2] /2 
s2 = [ (el+e4) /2-e2] /2 
s3 = -[e3- (el-e4) /2] /2 
s4 = [e3+ (el-e4) /2] /2 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = (el+e3)/2 
s2 = (e2-e4) /2 

s3 = [(el-e3).A - (e2+e4).B]/2 
s4 = [-(el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [(el-e3).A - (e2+e4).B]/2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
s3 = (el+e3) /2 
s4 = (e2-e4)/2 
This design is inferred from the design of Fig. 13 
by inverting values el, e2, e3, e4 and values, si, s2. 
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s3, s4, and by replacing B with -B because coefficient 
is now of the type e^*^"^^*^*. As A and B represent the 
cosine and sine of the same number, they have become 
A^+B^=l . The expressions of si, s2, s3, s4 are thus 
5 simplified. 

According to the invention, addresses are 
furthermore generated for sample addressing. Provision 
is made for generating two binary addresses of li bits 
per butterfly, each binary address corresponding to the 

10 rank of a butterfly input sample. The addresses of the 
other two samples to be applied to the butterfly inputs 
are obtained by inverting the p-p low-order bits of the 
first two addresses* 

In the same way as for the transformation methods 

15 operating according to time interleaving, the two 
binary addresses produced are consecutive for an 
internal butterfly. For a peripheral butterfly, the p-p 
low-order bits of the first generated address are equal 
to 0, and the yi-p low-order bits of the second address 

20 form a number equal to N/2^*^-l. By way of example, the 
addresses produced for the transformation circuit of 
Fig. 22 are gathered in Fig. 24. 

Advantageously, provision can be made for storing 
even address samples and odd address samples in two 

25 separate memories in order to reduce processing time of 
the transformation operation. 

Finally, the addresses 'generated for addressing 
the samples are also used for addressing coefficients 
W^. The value o'f the parameter s is used for addressing 

30 the corresponding coefficient W^. In this embodiment, 
the parameter s is equal to: 
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if p-p-l=p-2, the. number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, 

if vi-p-l<vi-2, the number formed by the ia-p-1 
5 low-order bits of the second address produced for said 
internal butterfly, followed by p-1. zero bits at the 
end of the number, 

if p-p-l>ii-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
10 internal butterfly, minus its p+1 low-order bits. This 
case corresponds to the first step (p=0) of the 
transformation methods operating according to a 
frequency interleaving algorithm- 




CIAIMS 

1, A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
5 digital signal' defined by a series of real starting 
samples x(n), with N a power of two and n e [0,.N-1], 
comprising supcessive transformation steps (2) for 
transforming input samples into output samples, all the 
transformation steps being performed by means of a 

10 single set of butterflies with several inputs and 
several outputs, the operating mode of which is 
modified selectively in each transformation step, the 
input and output samples of each transformation step 
being stored in a storage memory, a series of N output 

15 samples y (n) representative of the fast Fourier 
transform or the inverse fast Fourier transform of the 
output samples x(n) being provided in the last 
transformation step, 

characterized in that output samples y(h) are 

20 real, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x(n) processed in the first transformation step 
25 are classified in bit-reversed order of their index n, 
output samples y{n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 
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y(0) = Re[X(0)] 

y(n) = Re tX ( (n+1) /2) ] for n being odd and 

different from N-1 

Yin) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = Re[X(N/2) ] 
where samples X(n), with n e [0,.N-1], designate the 
complex samples of the series corresponding to the fast 
or inverse fast Fourier transform of the starting 
sample series x(n) , 

2. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N complex samples 
X(n) conjugated by pairs represented by a series of N 
real starting samples y(n)/ with N power of two and n 
e [O.-N-l], the starting samples y(n) being defined as 
follows: 

y(0) = Re[X(0) ] 

y(n) = Re[X ( (n+1) /2) ] for n being odd and 

different from N-1 

y{n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = RetX (N/2) ] 
this calculation method comprising successive 
transformation steps for transforming input samples 
into output samples, a series of N output samples x{n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 




is modified selectively in each transf orraation step, 
the input and output samples of each transformation 
step being stored in a storage memory, 

characterized in that output samples x(n) are 

5 real, 

and in that -the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
10 are classified in ascending order pf index n, the 
output samples x(n) are provided in the last 
transformation step in bit-reversed order of index n. 

3. The calculation method according to claim 1 or 
2, characterized in that, in each transformation step, 

15 each butterfly transforms input sample pairs, the ranks 
of the input samples of the same pair within the series 
of input samples of said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 

20, butterfly. 

4. The calculation method according to claim 3, 
characterized in that it comprises yi-l transformation 
steps Ep with ti=log2 (N) and p e [0..\i-2]. 

5. The calculation method according to claim 4, in 
25 turn dependent on claim 3, in turn dependent on claim 

1, characterized in further comprising: 

a preliminary step of modifying the sequence of 
the starting samples x(n) ranked in ascending order of 
index n and showing them in bit-reversed order of index 
30 n in the first transformation step, and 
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- a final step of processing the series of output 
samples y(n) and providing a series of N complex 
conjugated samples X(n) corresponding to the fast or 
the inverse fast Fourier transform of the series of 
5 starting samples x{n). 

6. The calculation method ,of claim 4/ in turn 
dependent on claim 3, in turn dependent on claim 1, or 
according to claim 5, characterized in that, in each 
transformation step Ep, butterflies are distributed 

10 among calculation blocks, 

in that each calculation block has a peripheral 
butterfly and/or 2^-1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
15 the input samples of rank 2^^^ar 2P^^a+2P''^-l, 2P^2a+2P^S 
2*^"^ 2^*^-1 into output samples of the same rank, 

and in that the internal rank t butterfly of the 
rank a calculation block in transformation step Ep 
transforms the input samples o-f rank 2'^'^^a+2x+l, 
20 2P^2^+2t+2, 2P"2a+2P^^-2x-3, 2P^2a+2P^2-2x-2 into output 
samples of the same rank, with p>l . 

7. The calculation method according to claim 6, 
characterized in that each butterfly is assigned a 
coefficient W^, whereon the calculation inside the 

25 butterfly is based, said coefficient being equal to e" 
j(2«s/N) ^^^Yl s e [0.. N/4-1] for a fast Fourier transform 
and is: equal to e^*^''^^'** with s e [0., N/4-1] for an 
inverse fast Fourier transform. 

8. Calculation method according to claim 7, 
30 characterized in that the internal rank x butterfly of 



the rank a calculation block in transformation step Ep 
is assigned coefficient with 5 = (x+l) . . 

9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 

5 the transformation steps are all of the same type and 
have 

four inputs for receiving input samples and 
four outputs ^or providing output samples, 

four additional inputs, respectively primary 
10 mode> secondary mode, permutation, and coefficient 
inputs, 

in order to selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
15 the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs, 

10. The calculation method according to claim 9, 
characterized in that, for each butterfly, the primary 

20 mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the other ones . 

25 11. The calculation method according to claim 10, 

characterized in that, in transformation step Ep, each 
calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies. 

12. The calculation method according to claim 11, 

30 characterized in that the secondary mode signal is 1 if 
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the peripheral butterfly is used for the last 
transformation step, and otherwise 0. 

13. The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
5 e3, and e4, and for a complex coefficient W^=A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

.1) if the primary mode and secondary mode signals 
are 0: si = el + e2 
10 s2 = el - e2 

s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

15 si = el + e2 + e3 + e4 

s2 = el " e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode signal is 1 and the 
20 permutation signal is 0: 

si = el + A.e3 - B.e4 

s2 = e2 + B,e3 + A.e4 

s3 = el - A,e3 + B.e4 

s4 = -e2 + B.e3 + A.e4 
25 4) if the primary mode signal is 1 and the 

permutation signal is 1: 

si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 

s3 = el - A. e3 - B • e4 
30 s4 = e2 + B.e3 + A,e4 



14, The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 
calculation block comprises: 

2^-1 internal butterflies and a peripheral 
5 butterfly for the even values of index p as well as for 
the last transformation step if ^ is even, and 
- 2^-1 internal butterflies, otherwise. 
_ 15. The calculation method according to claim 13, 
characterized in that the secondary mode signal is 1 if 
10 the peripheral butterfly is used for the last 
transformation step with \i being odd, and otherwise 0, 

16, The calculation method according to claim 15, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j .B, the 
15 butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
20 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1 : 

25 si = el + e4 

s2 = e2 
s3 = e3 
s4 = el - e4 

3) if the primary, mode signal is 0 and the 
30 permutation signal is 1: 

si - (e3 + e4) - (el + e2) 



s2 = el " e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A, e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B.e4 
s4 = -e2 + B.e3 + A.e4 

5) if the primary mode signal, is 1 and the 
permutation signal is 1 : 

si = el - A.e3 + B,e4 
s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 - B.e4 
s4 = e2 + B.e3 + A.e4 

17. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 
calculation block comprises: 

2^-1 internal butterflies and a peripheral 
butterfly for the even values of index p, and 
2^-1 internal butterflies, otherwise. 

18. The calculation method according to claim 17, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the first 
transformation step with p being even, and otherwise 0. 

19. The calculation method according to claim 18, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A-hj .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 
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1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 
s2 = el - e2 
5 s3 = e4 - e3 

_ - - ' -s4 = e3 e4 " ^ " - - -^-^ 

2) if primary mode, secondary mode and permutation 
signals are 0: 

si = el + e2 + e3 + e4 
10 s2 = el - e2 

s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
are 0 and the permutation signal is 1: 

15 si = (e3 + e4) - (el + e2) 

s2 el - e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 

4) if the primary mode signal is 1 and the 
20 permutation signal is 0: 

si = el + A.e3 - B.e4 

s2 = e2 + B.e3 + A.e4 

s3 == el ~ A.e3 + B.e4 

s4 = -e2 + B.e3 + A.e4 
25 5) if the primary mode signal is 1 and the 

permutation signal is 1: 

si = el - A.e3 + B.e4 
^ s2 = -e2 + B.e3 + A,e4 

s3 = el + A.e3 - B.e4 
30 s4 = e2 + B.e3 + A.e4 
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20 • The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - four inputs for receiving input samples and 

four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

10 in order to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 

15 additional inputs, 

and in that the final step furthermore performs an 
addition and subtraction between the first and the last 
output sample provided in the last transformation step. 

21. The calculation method according to claim 20, 
20 characterized in that, in transformation step Ep, each 

calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies. 

22. The calculation method according to claim 21, 
characterized in that, for four input samples el, e2, 

25 e3, and e4, and for a complex coefficient W^=A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
si = el + e2 
30 s2 - el " e2 




s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
5 s2 = e2 + B.e3 + A.e4 

s3 = el - A*e3 + B.e4 
s4 - -e2 + B.e3 + A,e4 
-3) if the primary signal is 1 and the permutation 
signal is 1: 
10 si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
s3 = el + A.e3 "B\e4 
s4 = e2 + B.e3 + A.e4 
23. The calculation method according to claim 9 or 
15 20, characterized in that the first and second binary 
addresses of |i bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 
20 24. The calculation method according to claim 23, 

characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

25. The calculation method according to claim 23 
25 or 24, characterized in that, if a peripheral butterfly 

is involved, the p+2 low-order bits of the first 
address are equal to 0, and the p+2 low-order bits of 
the second address form a number equal to 2^"'^-l. 

26. The calculation method according to claim 24 
30 or 25, characterized in that the address of the. two 

other samples to be applied to the inputs of the 



butterfly, be they peripheral or internal, are obtained 
by inverting the (p+2) low-order bits of said first and 
second produced addresses. 

27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories. 

-28. The calculation method according to claim 25, 
characterized in that the value of the parameter s of 
the coefficient assigned to an internal butterfly in 
transformation step Ep is coded by p-2 bits, and is: 

if p+l=p-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, 

if p+l<vi-2, the nximber formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, followed by p-p-3 zero 
bits at the end of the number, 

if p+l>vi-2, the number formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, minus its p-p-l low-order 
bits. 

29. The calculation method according to claim 4, 
in turn dependent on claim 3, in turn dependent on 
claim 2, characterized in that in each transformation 
step Ep, the butterflies are distributed among 2^ 
calculation blocks, 

in that each calculation block comprises one 
peripheral butterfly and internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 



the input samples of rank 2^'^a+2^"P'^-l, 2^'^a'^2^'^' 

^, 2^"'*a+2^~P-l into output samples of the same rank, 

and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
5 transforms the input samples of rank 2^'Pa+2x+l, 

Pa+2x+2, 2^-Pa+2^-P-^-2x-3, 2^-Pa+2^^-P-2x-2 into output 
samples of the same rank. 

30. The calculation method according to claim 29, 
characterized in further comprising a final step of 

10 modifying the sequence of the output samples provided 
in the last transformation step and classifying them in 
ascending order of index n. 

31. The calculation method according to claim 29 
or 30, characterized in that each butterfly is assigned 

15 a coefficient W^, whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j(2«s/N} ^^^^ 3 ^ [0.. N/4-1] for a fast Fourier transform 
and is equal to e^*^"^^^* with s € [0., N/4-1] for an 
inverse fast Fourier transform. 

20 32. Calculation method according to claim 31, 

characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 
is assigned coefficient with 5 = (x+1) .2^. 

33. The calculation method according to claim 32, 

25 characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

four inputs for receiving input samples and 
four outputs for providing output samples. 



four additional inputs, respectively primary 
mode/ secondary mode, permutation, and coefficient 
inputs, 

in order to selectively apply different 
5 transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 
additional inputs . 
10 34.. The calculation method according to claim 33, 

characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

in that the permutation signal is 0 for the even 
15 rank calculation blocks, including rank 0, and 1 for 
the odd values. 

35. The calculation method according to claim 31 
or 34, characterized in that the secondary mode signal 
is 1 if the butterfly, be it peripheral or internal, is 

20 used for the first transformation step, and otherwise 
0. 

36. The calculation method according to claim 35, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j .B, the 

25 butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: 

si = (el + e2)/2 
30 s2 = (el - e2) /2 

s3 = (e4 - e3) /2 



# so 0 



s4 = (e3 + e4)/2 
2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = [ (el+e4) /2-e2] /2 

s2 = [ (el + e4) /2-e2] /2 

s3 = [e3^(el-e4)/2]/2 

s4 = [e3+ (el+e4) /2] /2 
- 3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = (el+e3)/2 

s2 = (e2+e4)/2 

S3 = [(el-e3).A- (e2+e4).B]/2 
s4 = [-{el-e3).B + (e2+e4).A]/2 
4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [ (el-e3) .A - (e2+e4).B]/2 
s2 - [-{el-e3).B + {e2+e4).A]/2 
s3 = <el+e3)/2 
s4 = (e2-e4)/2 

37. The calculation method according to claim 33, 
characterized in that the first and second binary 
addresses of ]x bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 

38. The calculation method according to claim 31, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

39. The calculation method according to claim 37 
or 38, characterized in that, if a peripheral butterfly 



is involved, the vi"P low-order bits of the first 
address are equal to 0, and the p-p low-order bits of 
the second address form a number equal to 1^/2^*^-1. 

40. The calculation . method according to claim 38 
5 or 39, characterized in that the address of the two 

other samples to be applied to the inputs of the 
butterfly are obtained by inverting the p-p low-order 
bits of both produced addresses, 

41. The calculation method according to claim 40, 
10 characterized in that even-numbered address samples and 

odd-numbered address samples are stored in two separate 
memories. 

42. The calculation method according to claim 41, 
characterized in that the value of the parameter s of 

15 the coefficient assigned to an internal butterfly in 

transformation step Ep is coded by p-2 bits, and is: 

if ii-p-l=u-2, the number formed by the p-p-1 

low-order bits of the second address produced for said 

internal butterfly, 
20 - if ia-p-l<ia-2, the number formed by the p-p-1 

low-order bits of the second address produced for said 

internal butterfly, followed by p-1 zero bits at the 

end of the number, 

if ia-p-l>vi-2, the number formed by the p-p-l 
25 low-order bits of the second address produced for said 

internal butterfly, minus its p+1 low-order bits. 




ABSTRACT 

METHOD OF CALCULATING THE FAST FOURIER TRANSFORM AND 
THE INVERSE FAST FOURIER TRANSFORM 

The invention relates to a method of calculating 
the fast Fourier transform or the inverse fast Fourier 
transform of a series of N real samples x(n), with N 
power of two, operating according to a time 

10 interleaving algorithm and providing the sample series 
X(n) in ascending order of index n and using limited 
calculation and storage means. The invention also 
relates to a method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 

15 series of N conjugated complex samples X{n), with N 
power of two/ operating according to a frequency- 
interleaving algorithm. 
Fig. 8 

Application: image or acoustic signal processing, 
20 multicarrier modulation. 



2/17 






J 




3/17 





M 
Ul 




e4 
Ul 



ffi 
R 

s 

O 



o 



S2. 

< 



>3 




T 

























5/17 



• 



[X(0) ,X(1) ,X(2) ,X(3) ,X(4) ,X(5) ,X(6) ,x(7) , 
X(8) ,X(9) ,X(10) ,X(11) ,X(12) ,X(13) ,x(14) ,x(15) ] 



PRELIMINARY STEP 4 



Y 



[XCO) ,X(8) ,x(4) ,X(12) ,X(2) ,X(10) ,X(6) ,X(14) , 
X(l) ,X(9) ,X(5) >X(13) ,X(3) ,X(11) ,X(7) ,X(15) ] 







TRANSFORMATION. STEPS 


> 


f 



Cy(0) ,y(l) ,y(2) ,y(3) ,y(4) ,y(5) ,y(6) ,y(7) , 
y(8) /y(9) ,y(io) ,y(ii) ,y(i2) ,y(i3) ,y(i4) ,y(i5) ] 



i 



FINAL STEP 



i 



[X(0) ,X(1) ,X(2) ,X(3) ,X(4) ,X(5) ,X(6) ,X(7) , 
X(8) ,X(9) ,X(10) ,X(11) ,X{12) ,X(13) ,X(14) ,X(15) ] 



FIG. 6 



x(2)^-x(10)=A^(4) 



|x(6)>x(14)=/W6r 
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A«(4)+A«(6)=Br(4) 



iraj A«{4)-A«(6)=Br(6) 
I 




»H •[A«(4)-A<,(6)I=.Br(S) 



A«(4)>A^(6)=Baf4^ 



FIG.7A 




6«(9)-»-3.B,t(t 3)^3.St(1 3)=:Cr(9) 
B|(9)-3.B«0 3)+3.8»(l 3)=C{9) 



BR(9)-3.Bft(1 3)-3.B,(1 3)=Cr(1 1 ) 
-B,(9)-3.Bft(1 3)»3.Bi(13)=C,(t 1 } 




S«(9).3.B.^(1 3>-3.B,( 1 3)=Cn( 1 1 ) 



■B.(9)>3.8,i(13)»3.B.(t3)=C,(t t) 



8ft(9)»3.B^(t3)>3.Bt(13)=Cn(9) | 
8,(9)-3.B.,(1 3)>3.Bi(1 3)=Ct(9) | 
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FIG. 18 
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(2nd TRANSFORMATION STEP) 
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